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Abstract 

The effects of the finite charm quark mass on bottom quark mass determinations from T sum 
rules are examined in detail. The charm quark mass effects are calculated at next-to-next-to- 
leading order in the non-relativistic power counting for the T sum rules and at order for the 
determination of the bottom MS mass. For the bottom IS mass, which is extracted from the 
T sum rules directly, we obtain M^^ — 4.69 ± 0.03 GeV with negligible correlation to the value 
of the strong coupling. For the bottom MS mass we obtain Mb(Afb) = 4.17 ± 0.05 GeV taking 
ai"'"^'' (Afz) = 0.118 ± 0.003 as an input. Compared with an analysis where all quarks lighter 
than the bottom are treated as massless, we find that the finite charm mass shifts the bottom IS 
mass, Ml^, by about -20 MeV and the MS mass, Mb(Mb), by -30 to -35 MeV. 
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1 Introduction 



Precise and accurate determinations of the bottom quark mass parameter are mainly motivated by 
current and future B physics experiments, which aim at improved measurements of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements. This will hopefully get us closer to an understanding 
of the origin of flavour mixing, mass generation and CP violation. A precise and accurate knowledge 
of the bottom quark mass (in the scheme that is appropriate for the particular application) will 
eliminate one source of uncertainties. In particular, for the determination of Vub, and also Vet, from 
the semileptonic B decays the bottom quark mass is important in view of the strong dependence of 
the total rate and certain distributions on the bottom mass parameter. Another field where a precise 
bottom quark mass is desirable is represented by grand unification models. Here, precise values of the 
MS running heavy quark masses that are determined at low energies and that are as much as possible 
free of uncertainties from the strong interaction at low scales can constrain the allowed parameter 
space for a given scenario of flavour generation. 

The bottom quark mass parameter can be obtained from properties of hadrons containing 
one or more bottom quarks using perturbative quantum chromodynamics (QCD) or lattice methods. 
Recently, lattice calculations have been applied to bottom quark mass determinations from B and 
T meson spectra whereas perturbative methods using non-relativistic effective theories were 

applied to bottom quark mass determinations from the mass of the T(1S) and sum rules for the 
T masses and electronic widths 0, ^, ^, ^, ^, ^, |lO| . Extractions of the bottom quark MS mass at the 
scale Mz have also been carried out from analyses of LEP data on bottom-antibottom-gluon three-jet 
events [11|. The results of the bottom MS mass from these analyses have rather large uncertainties, 
but they have established the evolution of the MS bottom quark mass at higher scales. 

Using perturbative methods, non-relativistic sum rules for the masses and electronic decay 
widths of T mesons, bottom-antibottom quark bound states that have photonic quantum numbers 
and that can be produced in e^e~ annihilations, are in principle the most reliable tool to determine 
the bottom quark mass. Using causality and global duality arguments, one can relate energy integrals 
(or "moments") over the total cross section for the production of final states containing a bottom- 
antibottom quark pair in e^e~ collisions to derivatives of the vacuum polarization function of bottom 
quark currents at zero-momentum transfer. For a particular range of numbers of derivatives the 
moments are saturated by the very precise experimental data on the T mesons and, at the same time, 
can be calculated reliably using perturbative QCD in the non-relativistic expansion. In contrast to the 
bottom quark mass determination from the T(1S) mass, the sum rule analysis is practically free from 
model assumptions on non-perturbative effects as long as the number of derivatives that is employed is 
not too high. Because the moments have, for dimensional reasons, a strong dependence on the bottom 
quark mass, these sum rules can be used to determine the bottom quark mass to high precision. 

In the most recent sum rule analyses at NNLO in the non-relativistic expansion ||8|, ^, IC] 
special care was taken to consistently eliminate the strong linear sensitivity to small momenta and 
the associated (artificially) large perturbative corrections that were contained in previous analyses 
that employed the bottom quark pole mass. This was achieved by implementing bottom quark short- 
distance mass definitions into the moments that were constructed, particularly for the situation where 
the bottom quark is very close to its mass shell. For the implementation of these masses the systematic 
cancellation of the linear infrared-sensitive contributions and the compliance with the non-relativistic 
power counting had to be ensured in each order of perturbation theory. Melnikov-Yelkhovsky Q used 
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the kinetic mass, Hoang |p the IS mass and Beneke-Signer |1C] the PS mass. These "low-virtuahty 
short-distance masses" can be used directly for other processes where the bottom quark virtuahty is 
small with respect to its mass. In semileptonic B meson partial rates, for example, they lead to a 
very well behaved perturbative expansion (see e.g. Ref. ||9|, |l^]). For situations where the bottom 
quark is very off-shell, on the other hand, the MS mass definition is the appropriate one. Starting 
from the low-virtuality short-distance bottom masses the previously mentioned analyses obtained 
Mb(Mb) = 4.2 ± 0.1 GeV §, 4.20 ± 0.06 GeV |] and 4.25 ± 0.08 GeV @ for the bottom MS mass. 
These results represent determinations of the bottom MS mass at order in the usual loop expansion 
of perturbative QCD. The results are compatible with each other and also with the most recent lattice 
calculations (see Ref. |Q] for a convenient compilation of results obtained after 1994). The uncertainties 
are considerably smaller than Aqcd and amount to about 2%. It will be the aim of future analyses to 
further reduce the uncertainty, to maybe 1% or less, i.e. to an amount of a few times 10 MeV. Thus, 
it is mandatory to account for effects that can lead to for shifts of several times 10 MeV. 

One effect, which has been entirely neglected in all previous sum rule analyses, is coming from 
the masses of those quarks that are lighter than the bottom, most importantly the charm quark. The 
massless approximation for the light quarks is a good one, provided that the light quark mass is much 
smaller than any of the relevant dynamical scales. In this case it is legitimate to simply expand in 
the light quark mass. Because all linear sensitivity to small momenta is eliminated by the use of a 
short-distance bottom quark mass, the first non-vanishing correction to the massless limit is of order 
(as/7r)^m^jgjj^./Mb, where mught is the light quark mass. Even for the charm quark this would amount 
to a shift of only a few MeV in the bottom quark mass and could be safely neglected. However, 
in order to judge whether the massless approximation is really justified, one needs to compare the 
light quark masses with the scales that are relevant to the non-relativistic bottom-antibottom quark 
dynamics that is encoded in the moments. There are two issues that make the effects of light quark 
masses, and of the charm quark mass in particular, a subtle problem: one is specific to the extraction 
of the bottom quark MS mass from the non-relativistic sum rules and the other exists for any process 
that is calculated at large orders of perturbation theory. 

The issue that is specific to the dynamics of the non-relativistic bottom-antibottom pair en- 
coded in the moments arises from the fact that the inverse Bohr radius, M^,as, is a relevant dynamical 
scale. Because the inverse Bohr radius is about the same as the mass of the charm quark, it is a priori 
not allowed to expand in the charm quark mass. Rather, the effects of the charm quark mass have 
to be treated exactly, leading to a non-trivial dependence on the ratio of the charm quark mass and 
the inverse Bohr radius. In Ref. [13| is has been shown that this leads to an incomplete cancellation 
of linear charm quark mass terms oc as^richarm in the relation between the bottom MS and IS masses. 
The resulting shift in the MS mass Mb(Mb) was estimated to be around —20 MeV, based on a NLO 
calculation (which corresponds to an order accuracy in the MS mass). If an expansion in the charm 
quark mass would be allowed, all linear charm quark mass terms would cancel exactly. The large size 
of the shift caused by the incomplete cancellation of the linear charm quark mass term comes from 
the fact that the linear charm quark mass term is enhanced. Technically, this happens because linear 
light quark mass terms represent non-analytic terms (in the square of the mass) that are multiplied by 
a factor tt^, and because the scale of the strong coupling governing them is the light quark mass and 
not the bottom mass. Physically, this happens because light quark masses act as an infrared cutoff for 
the gluon line in which the light quark loops are inserted. A linear sensitivity to small momenta is, 
after expansion, reflected by a non-analytic, linear dependence on the light quark mass. Because this 
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linear, non-analytic dependence can only arise from momenta of the order of the light quark mass, 
the strong coupling multiplying the linear terms is renormalized at the light quark mass. Thus, the 
masses of the light quarks work in complete analogy to the infinitesimally small fictitious gluon mass 
that is sometimes used in standard renormalon analyses. The difference is that the light quark masses 
are not fictitious, but finite numbers that are given by experimental data. Whether we are allowed to 
expand in them (or not) depends on their relation to the dynamical scales of the problem. 

The other, more general issue is associated with the fact that QCD perturbation theory at 
large orders of perturbation theory becomes dominated by exponentially decreasing momenta. This 
is the origin of the divergent asymptotic behaviour of QCD perturbation theory associated to the 
so-called "infrared renormalons" |14, 15, 16 1. This means that the genuine effects from light quarks 
with finite masses effectively decouple, as from a certain order, and that the large order asymptotic 
behaviour is considerably affected by the masses of the light quarks 0]. In fact, at large orders 
of perturbation theory, the difference between the massless approximation for light quarks and a 
calculation that includes the light quark masses becomes arbitrarily large because the evolution of the 
strong coupling differs in the two cases for momenta of order or below a specific light quark mass. 
Even at intermediate (i.e. phenomenologically relevant and calculable) orders, however, the effective 
decoupling of light quarks can lead to considerable deviations from the massless approximation, if the 
perturbative coefficients are already dominated by sufficiently small momenta. The effective decoupling 
is therefore not just a purely academic issue. 

From the issues mentioned above we can expect the following qualitative effects coming from 
the finite charm quark mass in the bottom MS mass obtained from the bottom IS mass: the largest 
correction arises at order from the fact that one can expand in the charm quark mass in the 
bottom pole-MS mass relation, but not in the pole-lS mass relation |l^. The higher order charm 
quark mass corrections should decrease because the importance of the dynamical scales M^, and 
Mbtts becomes less prominent as the dominating momenta decrease exponentially. This means that 
the charm mass corrections are under control and can be calculated reliably^]. From the light quarks 
other than the charm, we do not expect any sizeable corrections because, at the accessible orders of 
perturbation theory, their masses are simply too small to make any of the previous considerations 
relevant. However, a quantitative statement about the actual size and behaviour of the charm mass 
corrections can only be made through an actual calculation. In addition, is it mandatory to assess 
how the finite charm quark mass affects the extraction of the bottom IS mass from the T sum rules. 
For the bottom MS mass the corrections might either add up or cancel at the end. In the former case 
the charm quark mass effects for the final value of Mb(Mi,) might be larger than the 20 MeV estimate 
that has been obtained in the NLO (or order a^) analysis of Ref. [13|. 

In this paper we determine the charm mass corrections at order in the bottom MS-IS mass 
relation as well as for the bottom IS mass determination from the T meson sum rules at NNLO in the 
non-relativistic expansion. This includes the calculation of the NNLO light quark mass corrections 
to the perturbative contributions of the T(1S) mass in terms of the pole mass, and the order 
light quark mass corrections in the heavy quark pole-MS mass relation. The latter corrections are 
determined in the linear light quark mass approximation because an expansion in the light quark 
masses is allowed in that case. The determination of the linear order light quark mass corrections 

^ We emphasize that this statement is not correct if one attempts a determination of the bottom quark pole mass. Its 
weU-known ambiguity of order Aqcd is directly reflected in an uncontrollable behaviour of the charm mass corrections. 
This is demonstrated in Sec. 5.2. 
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is based on the conjecture that linear hght quark mass terms are absent in the total static energy 
of a heavy quark-antiquark pair. As far as the bottom MS-IS mass relation is concerned, we also 
attempt to estimate the order charm mass corrections based on a calculation of the NNNLO 
(N^LO) perturbative corrections in the T(1S) mass, in the large-/3o approximation. For the T sum 
rules we determine charm mass corrections at NLO and NNLO in the non-relativistic expansion, where 
we neglect the double-insertion contributions coming from second order Rayleigh-Schrodinger time- 
independent perturbation theory at NNLO. For all results we discuss in detail the behaviour of the 
light quark mass corrections. It is demonstrated that the effects coming from the masses of up, down 
and strange quarks expected, be neglected. 

We note that, although the effects of the light quark masses can be presented in a straight- 
forward way in the framework of effective theories for non-relativistic heavy quark-antiquark pairs, 
we will not give such a presentation in this work. This is because we feel the discussion of the many 
interesting aspects of light quark mass effects, which do not depend at all on the concepts of effective 
theories, would become cumbersome. We also emphasize that all results presented in this work are 
given in the convention that the evolution of the strong MS coupling includes the light quarks with 
finite mass, i.e. we use ai""'"^^ for our results for the bottom quark masses. We do this because all 
previous analyses, which neglected light quark masses, naturally used this definition of as- It is there- 
fore easy to update the previous analyses by simply adding the new corrections determined in this 
work. In this context we also note that many conclusions and statements made in this work about the 
behaviour and the size of the charm mass corrections at the different orders are not adequate if ai"'~'^^ 
is chosen as the definition of the strong coupling, since this would lead to a reshuffling of corrections. 
The final numbers for the bottom IS and MS masses, however, do not depend on this choice. 

The plan of this paper is as follows: in Sec. |2| we briefiy review the ingredients needed to 
describe the heavy-quark-antiquark dynamics at NNLO in the non-relativistic expansion for massless 
light quarks and the result for the NNLO light quark mass corrections to the static potential obtained 
earlier by Melles 19 1. The results by Melles are rewritten in a more convenient dispersion relation 
representation. The limit of vanishing light quark mass and the resulting linear light quark mass terms 
are examined, and the prediction for the order linear light quark mass corrections in the heavy 
quark pole-MS mass relation is made. The quality of the linear light quark mass approximation 
in the heavy quark pole-MS mass relation is analysed in Sec. ^. In Sec. |^ the NNLO light quark 
mass corrections to the relation between the heavy quark IS and pole masses are presented, and 
in Sec. ^ the results from Sees. |3| and ^ are combined to derive the order (or NNLO) relation 
between the heavy quark MS and IS masses, where light quark mass effects are taken into account. 
The upsilon expansion a prescription for a modified perturbative expansion that allows for a 
consistent derivation of the MS-IS mass relation in each order, is reviewed. The size of the charm 
mass corrections is analysed and an estimate for the order (or N'^LO) charm mass corrections 
is given. In Sec. ^ we determine the NNLO light quark mass corrections to the moments of the T 
sum rules, and in Sec. |^ we carry out the numerical analysis to determine the bottom IS mass from 
fitting theoretical and experimental moments. The relation of our result for the bottom IS mass to 
the mass of the T(1S) meson is discussed, and the result for the bottom MS mass is presented. All 
results are discussed focusing on the charm mass effects. Section ^ contains a summary. Attached 
to the paper are five appendices: in App. ^ details on the calculations of the NNLO light quark 
mass corrections to the mass of the J^'-' = 1 , ^Si heavy perturbative quarkonium ground state are 
given. Appendix ^ contains the calculation of the N'^LO perturbative corrections to the mass of the 
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heavy quarkonium ground state in the large-/3o approximation. In App. ^ the results for the order 
(or N'^LO) perturbative corrections to the heavy quark MS-IS mass relation are derived in the 
large-/3o approximation. Appendix ^ contains exact analytic results for the order light quark mass 
corrections to the heavy quark pole-MS mass relation coming from vacuum polarization insertions, 
and App. ^ gives some formulae for the corrections to the moments of the T sum rules from massless 
light quarks. 



2 Considerations on the Static Potential 
including Light Quark Mass Effects 

The static potential represents the central quantity needed to determine the effects of the masses of 
light quarks on the dynamics of a perturbative non-relativistic heavy-quark-antiquark pair. In this 
section we discuss the corrections to the static potential up to NNLO in the non-relativistic expansion 
coming from the mass of a light quark. 

2.1 Massless Light Quarks 

At NNLO in the non-relativistic expansion, the dynamics of a perturbative heavy-quark-antiquark 
pair^ represents a pure 2-body problem and can be described by a textbook-like Schrodinger equation 
that contains the kinetic energy term and a time-independent instantaneous potential. In configuration 
space representation, the Schrodinger equation has the form (E = y/s — 2Mq°'''): 



+ 



E 



X G{r,r\E) = 5^^\r - r') , (1) 



where Mq°'° is the heavy quark pole mass and s the squared centre-of-mass energy. Although the 
use of the pole mass is a bad choice for a phenomenological analysis, because of its strong sensitivity 
to small momenta, we will use it for intermediate calculations because it leads to the simple form of 
the Schrodinger equation shown in Eq. (|l|). The terms V^^'^ and V^^nlo j-gp^gggnt the static 

(Coulomb) potential at LO (order ~ a^), NLO (order v"^ ~ q^) and NNLO (order ~ a^), 
respectively, in the non-relativistic power counting, while lap is the Breit-Fermi potential known from 
positronium and Vna is a non-Abelian potential coming from a non-analytic gluonic correction to the 
exchange of a Coulomb gluon; Vbf and I^a are of NNLO. The kinetic energy terms — V^/Mq°'° and 
— V^/4(Mq°''')^ are of LO and NNLO, respectively. Treating all ni quarks that a lighter than the 
heavy quark as massless the individual potentials read (a^ = a^"'\^), r = \r\, jl = e^^ ^, Ca = 3, 
= 4/3, r = 1/2) 

V^ir) = (2) 



An 



2 Po ln(/i r) -I- ai 



(3) 



^ In this context "perturbative" means that the heavy quark mass, Mq, is large enough, such that the hierarchy 
Mq 3> Mqv 3> Mqv^ S> Aqcd holds, V being the heavy quark velocity. 
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NNLO 
c.massless 



(r) 



(1^) ' [ /3o' ( 4 ln2(^ r) + y ) + 2 (2 /3o ai + /3i) ln(/i r) + 



(4) 



Cp UsTT 



1 + 3^^5-, 



5(3) (r) + 



Ci;' Us 1 
2(M^°"=)2 7 



+ (r V)V 



+ 



SCiT'Os 1 



2(M, 



Q 



bJ ' 



(5) 



V^NA(r-) 



(6) 



where 



11 4 



(7) 



Pi 



34 

y 



20 



Ca Tni — A:CfT rii , 



(8) 



31 ^ 20 



(9) 



a2 



4343 , 2 

h47r2 

162 



22 



- + C3 )Cl 



55 \ / 20 



1798 56 



(10) 



iS;, and are the bottom and antibottom quark spin operators, L is the angular momentum operator, 
and 7e = 0.57721566 ... is the Euler-Mascheroni constant. The constants ai and 02 have been 



calculated in Refs. |21, 22] and Refs. [23, 24|, respectively. For the determination of the corrections 
coming from a finite light quark mass in physical heavy-quark-antiquark properties at NNLO, the 
expressions for and K^^m^Ssioss '^i^l be needed explicitly. The superscript "(n;)" shown by the 
strong coupling indicates that the results are written in terms of the strong coupling that evolves with 
ni active quark flavours. If the light quark masses are much smaller than the typical inverse distances 
of the heavy quarks, this convention is the only sensible one, as it avoids the appearance of logarithms 
of the light quark masses. 



2.2 Effects of a Light Quark Mass 

From the explicit dependence of Eqs. ©-(P) on the number of light quarks n^, one can easily see that 
only the static potential at NLO and NNLO acquire additional corrections from a finite light quark 
mass. In the following we present the corrections to the static potential at NLO and NNLO for the 
case that only one light quark species has a finite mass m, whereas the ni — \ other quark species are 
treated as massless. For convenience we adopt the pole mass definition for the light quark mass m. 
The generalization to other light quark mass definitions, and to the case where more than one of the 
ni light quarks have a finite mass, is straightforward. We note that all final results for the light quark 
mass corrections presented in this work will be expressed in terms of the light quark MS mass m(m) . 
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Figure 1: NLO contribution to the static potential coming from the insertion of a one-loop vacuum 
polarization of a light quark with finite mass. 



At NLO the light quark mass corrections to the static potential arise from the insertion of the 
light quark self energy into the gluon line of the LO static potential (see Fig. 0). In momentum space 
representation the NLO static potential including the light quark mass corrections reads (a^ = ai"'^(/i)) 



NLO 



V 



(11) 



[12) 



where 



fix) 









(in 4 




) = 



















1 

3 ~ 
2p' 



{l-2z 

no 

dx f{x 



2 - Vl + 4z In 



/ V1 + 4Z + 1 
Vl + 4z- 1 



\fxP- 



11 + 



2x2 



(13) 
(14) 

(15) 



(16) 
(17) 



We note that P also contains the subtraction of the light quark mass vacuum polarization in the 
limit m ^ 0. Thus the light quark mass corrections 5V^^,j° vanish for m — > or ^ oq^ This is 
a consequence of our choice to include the light quark with mass m into the evolution of the strong 
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coupling. The dispersion relation representation for the light quark vacuum polarization function 



shown in Eq. (16) is quite useful because it allows for a quick determination of the light quark mass 



corrections in configuration space representation (m 



=7e 



m): 



cm V / 



V} 



dxf{x) 



~2mrx 



+ I ln(m r) + 



(18) 



The expression in the second line of Eq. (|lq) (without the subtraction of the contributions for m ^ 0) 



has been derived a long time ago by Serber and Ueling |25] for the effects of the electron on the 
QED static potential of a charged particle. The use of the dispersion relation representation for the 
light quark mass corrections to the static potential is advantageous for the determination of higher 
order terms in Rayleigh-Schrodinger perturbation theory due to the universality and simplicity of the 
dependence on r (or p). The remaining dispersion integration can be carried out numerically (or if 
possible analytically) at the very end. We believe that this method is actually the only feasible one 
to determine the light quark mass corrections, particularly for the light quark mass corrections in 
the NNLO static potential and for multiple insertions of the NLO static potential at higher order in 
Rayleigh-Schrodinger perturbation theory. For the determination of the dispersion relation represen- 
tation, it is a useful fact that the corrections of the light quark mass with the (pole) mass m decouple 
(i.e. vanish for p^ 0) if we exclude it from the evolution of the strong coupling. Up to two loops 
the matching relation of the strong coupling in the two schemes, where the light quark with mass m 
is included and excluded in the evolution, reads [E^] 



vr 

(ni-l) 



+ 



7 19 , / /i 

\ In — 

24 24 Vm 



2^ ^ 36 



56 V 



(19) 



Thus, after rewriting the static potential in terms of ai"' ^\ only the subtracted dispersion integral 
expression remains for the corrections coming from the finite light quark mass m. We have used this 
useful fact to determine the dispersion relation representation for the light quark mass corrections to 
the static potential at NNLO. 

At NNLO the light quark mass corrections to the static potential arise from dressing the one- 
loop Feynman diagram displayed in Fig. || with additional gluon, ghost and light quark lines, where 
the heavy quark lines in the resulting loops and the couplings of the heavy quarks to gluons correspond 
to those of static quarks [O] as they are used in Heavy-Quark Effective Theory. The corresponding 



diagrams were calculated numerically by Melles [18|. Fitted approximation formulae for the NNLO 
static potential with the corrections of a light quark mass in momentum as well as configuration space 
representation were given in Ref. |l^. There, the results were presented for n; = 1 and using ai"'~^^ 
as the definitions for the strong coupling. However, these results are somewhat cumbersome to use 
for bound state calculations at higher order in Rayleigh-Schrodinger perturbation theory. Using the 
decoupling properties of the massive light quark corrections as explained above, it is straightforward 
to derive the following dispersion relation representation of the light quark mass corrections to the 



10 



NNLO static potential in momentum space representation^ (a^ = ai"'^(/i)): 



where 



^NNLO / ) 

c,masslcss 



47r. 
a, \2 



TP 



m 



3 vp- 

4 x2 



/3oln(^) 



+ ai 



1 r 2 / 

/m \ \ ( 



In- 



161 
114 



26 
19 



(20) 
(21) 



(22) 



(23) 



dx 



1 



dx 



2cip' / 2 I .1 2 2 + ^ ^1 P / 2,^22 



4^2 



(24) 



and the definition of P has been given in Eq. (15). In contrast to the approximation formulae given 
in Ref. jl^, the contributions involving the one-loop vacuum polarization contributions coming from 
the light quark with mass m are exact. The remainder, which is parametrized in terms of the function 
X, depends on the constants ci, C2, di and d2- Because H vanishes for — > and X for p^ ^ 00, 
the constants satisfy the conditions ci + di = 1 and ci ln(4c2) + di ln(4(i2) = + ffCs- ^ useful 
parametrization for the constants, which fulfils these conditions and keeps C2 and d2 as free parameters. 



IS 



ci 



In A 
In^ ' 



In £2. 
In^ 



For the constants C2 and d2 we obtain 



/ 161 13 , , 
^ = ^^Pf 228 + 19 ^^-^"' 



C2 



do 



0.470 ±0.005, 



1.120 ±0.010 



(25) 



(26) 



(27) 



(28) 
(29) 



I thank M. Melles for checking Eqs. (23) and (E 
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from fitting to the results in Eq. (30) of Ref. |19]. We assume a 1% uncertainty for C2 and ^2, but 
we note that this estimate is not based on a statistical analysis. The uncertainties have been fixed 
by hand in an ad hoc manner taking into account the uncertainties quoted in Ref. and the fact 
that the numerical result shown in Eq. (30) of Ref. [|l9| ] does not decouple exactly. We note that the 
uncertainties in the charm mass effects for the bottom mass extractions caused by the errors in C2 and 
d2 are an order of magnitude smaller than those caused by the other approximations made in this work. 
We therefore ignore the errors in C2 and d2 for the rest of this paper. Finally, the NNLO light quark 
mass corrections to the static potential in configuration space representation read (a^ 

a, x2 
^37r^ 



<5K 



NNLO 
c,m 



(JXJ 

3 f 

- J dxf{x)e~ 



2mrx 



( A:vr?x^ 
Po ( m 7^ Ei(2 mrx) 



Ei(— 2 mr x)) — ai 



+ 3 In(mr) H — 
6 



/3oln(/ir) + ^j+/?o^ 



dx f{x) e 



- 2mrx 



5 1 

77 + ^ 



1 + — Vx"^ - 1 (1 + 2x2) In ( - 
2x V' 



+ Vx"^ - 1 
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where Ei is the exponential-integral function and T the incomplete gamma function, where r(0, ca) = 
joo ^-ax the exponential-integral function the relation Ei(z) + Ei(— z) = P /g°° 2t e^^'''*^^/(l — 

i^) dt^ where P refers to the principal value prescription, is quite useful. We note that of the three 
terms in the brackets on the RHS of Eq. (pO|), each one vanishes individually for m ^ 0. The first 
two terms incorporate all the light quark mass corrections that can be obtained from insertions of 
one-loop light quark vacuum polarizations. 



2.3 Small-Distance Limit and Linear Light Quark Mass Terms 

It is instructive to consider the light quark mass corrections to the static potential for the case where 
the distance between the heavy quark is much smaller then 1/m, i.e. in the limit m ^ 1/r. In this 
case the light quark mass corrections to the static potential read (a^ = ai"'^(/x)): 



m<^l/r 



ai + Po[ ln-^ + 3-41n2 



76 



301n2 ] + — [ciC2 + did2 



+ O(m^r) 



12 



4 /59 ^1 ^\ 76 , 



+ 0{m'^r), (31) 



where we have used the relation ai = |/3o — 8 in the second hne of Eq. (|3l|) . The order af/?o term 
in the second hne has aheady been derived in Ref. ||l^. The leading term in the expansion, which is 
linear in the light quark mass^, has a number of remarkable features: it is independent of r, enhanced 
by a factor of vr^ , and its characteristic scale is the light quark mass m rather than the inverse of the 
distance r. We emphasize that the linear light quark mass term cannot be obtained by first expanding 
Eqs. (|18|) and (^0|) in the light quark mass and carrying out the dispersion integrations (over the 
variable x) afterwards, because it contains contributions coming from large values of x, i.e. from 
momenta of the order of the light quark mass or smaller (see Eq. (p!^)). This means that the linear 
mass term is of infrared origin, which explains that it represents a non-analytic expression in m^. In 
fact, the linear light quark mass term is in complete analogy with the linear gluon mass term that 
would arise if we were give the gluon a infinitesimally small fictitious mass A. The linear light quark 
mass term at order a" directly corresponds to the linear gluon mass term that would be obtained at 
order a"^^. (This can be seen from Eq. ([l6|), using the change of variable = Am? x"^.) Just as for 
the linear gluon mass term, the existence of the linear light quark mass term in Eq. (^Tj) indicates 
that the static potential is, like the heavy quark pole mass, linearly sensitive to small momenta. Thus, 
the static potential also contains an (r-independent) ambiguity of order Aqcd [29|. This also explains 



why the perturbative series for the linear quark mass term is completely out of control for all light 
quarks, although the perturbative expansion of the static potential is most trustworthy in the small 
distance limit. For the charm quark the linear mass term amounts to 60 MeV at order and to 
86 MeV at order for mijn) = fj, = 1.5 GeV, = 0.36 and ni = A indicating a very bad higher 
order behaviour^. (See Sec. ^ for a more detailed discussion.) 

Fortunately, a linear sensitivity to small momenta and the corresponding ambiguity of order 



■Aqcd , is not contained in the total static energy |3^, |3^ 



Sstat = 2M^°'° + K(r). (32) 

In Refs. 32| this was proved at one loop. The proof in Ref. [31| demonstrates this also at the two- 
loop level, which is sufficient for the purpose of this work. A similar proof for inclusive semileptonic 



B meson decays can be found in Ref. |33]. Equation ( |3^ ) means that — if we are actually allowed 
to expand the static potential in the light quark mass — the linear light quark mass corrections are 
cancelled, if a short-distance mass definition is employed to describe the non-relativistic heavy-quark- 
antiquark dynamics. This remarkable feature allows us to determine the linear light quark mass 
corrections in the relation between the heavy quark pole and a short-distance mass, such as the MS 



* In a non-relativistic effective field theory for heavy quarks, in which the light quark is also integrated out, the linear 
light quark mass term corresponds to a residual mass term operator [^. This term accounts for the difference in the 
analyticity structure at low energies of the effective theory in which the light quark is integrated out and the theory in 
which the light quark is not integrated out. This also means that the heavy quark pole masses in both theories differ by 
the coefficient of the residual mass operator |2^], see also the section on quark masses in Ref. psf. 

® This statement about the higher order behaviour of the linear charm quark mass terms does not depend on which 
convention is used for the strong coupling. 
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one, up to order ap. 



1 

~ 2 



linear in m 



(33) 

linear in m 



The order linear light quark mass term of the static potential shown in Eq. ( |3T|) confirms the 
order hnear hght quark mass term in the heavy quark pole-MS mass relation calculated earlier in 
Ref. [^]. The order term that can be derived from Eqs. (|3l| ) and (33) has to our knowledge not 



been obtained before. A detailed discussion on the linear light quark mass corrections in the relation 
of the heavy quark pole and MS masses is given in Sec. ^. 

We are now in a position to substantiate the scenario for the light quark mass corrections that 
we have pictured in the introduction to this work for the extraction of bottom quark masses from T 
mesons at intermediate orders of perturbation theory. For the light quarks up, down and strange, the 
badly behaved linear mass terms indeed cancel in an extraction of the bottom MS mass, because their 
masses are much smaller than the dynamical scales Mb and (1/t) ~ MbQ^ and one can expand in them 
everywhere, in the same way as in the infinitesimally small gluon mass. For the charm quark, however, 
the expansion in its mass is allowed in the bottom quark pole-MS mass relation (because it is only 
governed by the scale Mq), but not in quantities involving the bottom-antibottom non-relativistic 
dynamics described through the Schrodinger equation (|^). This means that, for an extraction of 
the bottom MS mass from T mesons, at least some part of the large linear charm mass corrections 
remains uncancelled. Regarding the very bad perturbative convergence for the linear charm mass 
terms in Eq. ([3^), it is now natural to ask the question of how the remaining linear charm mass 
corrections will behave at higher orders. If they were to behave similarly to the linear charm mass 
corrections in Eq. (^T|) , then the precision of the bottom MS mass would be forever limited to several 
tens of MeV. Fortunately, this is not the case because the contributions from low momenta increase at 
higher orders of perturbation theory. This means that the small momentum region, where the linear 
mass terms originate, dominates over the larger momenta. Thus, although it is formally impossible 
to expand in the charm mass divided by the inverse Bohr radius, the coefficients of higher powers of 



this ratio will decrease at higher orders, making the linear term the dominant one. In Sec. 5.2 we will 
demonstrate this interesting feature explicitly. For the extraction of the bottom MS mass this means 
that the charm mass corrections are well under control and calculable. 



3 Light Quark Mass Corrections in the 
Heavy Quark Pole— MS Mass Relation 

The perturbative relation between the pole and the MS mass definition of a heavy quark is a main 
ingredient in the determination of the bottom MS mass from data on the T mesons. The mass of 
a light quark leads to corrections in this relation starting at order a^, originating from the insertion 
of a light quark vacuum polarization into the one-loop gluon line. For a proper determination of the 
bottom MS mass using calculations for the non-relativistic bottom-antibottom quark dynamics at 
NNLO (see Eq. (|l])), we need the bottom pole-MS mass relation at order af. Here, "proper" means 
that the artificially large perturbative contributions coming from a linear sensitivity to small momenta 
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are eliminated consistently in each ordeiQ For massless light quarks the heavy quark pole~MS mass 



relation is known analytically at order a^. The order corrections have been determined in Ref. |34] 
and the order corrections in Ref. |12]. A numerical determination of the order corrections can 



be found in Ref. |35|. In Ref. |34] the order corrections coming from the mass of a light quark have 
been determined analytically for all values of the light quark mass. The light quark mass corrections 
at order are only known in the linear mass approximation and have been determined in Sec. |2.3| of 
this work. At order the relation between the pole and the MS mass of a heavy quark for rii light 



quarks, of which one has a finite pole mass m, reads (a^ 



a 



M^{M^) 1 



5«( 



6^^^ ia 



(^q)): 

)+5iiliveKMQ(MQ),a, 



(3) TI (TI 
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where the corrections for massless light quarks read 

^(1) 



4 /Os 
3 V vr 



.(2) 



?(3) 



(^) ( 13.443 - 1.0414 r 
(^)^ ( 190.4 - 26.66 nz +0.6527 n 



(34) 

(35) 
(36) 
(37) 



The numbers displayed in Eq. (p7|) have been taken from Ref. |12|. The light quark mass corrections 
in the linear approximation read 
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If the MS definition is used for the charm quark mass, m = m(m)[l + |(^)], we have to apply the 
replacement 



f(2) 1^ TI (TJ 



-'massive 



m,M^{M^),as 



^(2) (m(m),MQ(MQ),a,) + e5«(a,)(5ili.,(m(m),MQ(MQ 



■'massive 



(2) t^(^\ TJ (TI 



(40) 



in Eq. (34). The generalization to the case where more than one of the light quarks has a finite mass 
is straightforward. The parameter e (e = 1) is an auxiliary variable that facilitates the determination 

^ The relation between perturbative orders (LO, NLO, etc.) in the non-relativistic power-counting scheme and 
perturbative orders {as, a^, etc., or one-loop, two- loop, etc.) in quantities unrelated to non-relativistic quark-antiquark 
physics, such as the pole-MS mass relation, can also be seen from Eqs. (|l|) and (||). The order linear light quark 
mass term is obtained from the static potential at NLO, the order linear light quark mass term is obtained from the 
static potential at NNLO, etc. 
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of the heavy quark MS-IS mass relation in Sec. |5|, where the non-relativistic power counting and the 
ordinary power counting via the number of loops have to be combined in such a way that the linear 
sensitive contributions cancel systematically in each order. Equation (^4|) is the basic equation that 
we will use for the determination of the bottom MS mass in this work. For the corrections coming 
from the charm quark mass, we will only use the linear mass approximation. We show in this section 
that this approximation should be sufficient to give a description of the charm mass effects with a 
relative uncertainty of 10%. We also demonstrate that the linear mass approximation becomes better 
at higher orders by explicitly comparing the linear mass terms with known exact results. 

At order the light quark mass corrections can be represented as an integral over the one-loop 
gluon virtuality times the subtracted light quark vacuum polarization P |3^, in analogy to the NLO 

r = m/Mq): 



light quark mass corrections to the static potential in Eq. (\\_ 
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where P is defined in Eq. (|T5| ) and 
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(43) 



Equation (|4l| ) is ultraviolet-finite owing to the subtraction of the massless vacuum polarization in the 
function P. From the first line of Eq. (|4l| ) it is straightforward to see that in the limit m <C Mq two 
momentum regions contribute to the integral: q ^ m and q ~ Mq. Thus, the expansion in terms 
of m/Mq cannot be obtained by naively Taylor-expanding the momentum integral in the first line of 
Eq. (|4l|), since this would only extract the contributions of the high momentum region q ~ Mq. The 
linear light quark mass term and all the other terms that are non-analytic in m?, such as and 
Inm, originate from the low momentum region q ^ m. To obtain the linear light quark mass 
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2n 



term directly, we therefore have to carry out a non-relativistic expansion in q^ /M^ 
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m/MQ 


0.1 


0.2 


0.3 


0.4 


0.5 


r(2),full /r4/a,N,2] 

^massive / L 3 V TT / -1 


0.117 


0.223 


0.320 


0.411 


0.496 


'-^massive/ L 3 ^ TT ' -1 


0.123 


0.247 


0.370 


0.493 


0.617 


^massive / ^massive 


1.06 


1.11 


1.16 


1.20 


1.24 



Table 1: Comparison of the full result (second line) and the linear mass approximation (third line) 
for the order a1 light quark mass corrections in the heavy quark pole-MS mass relation for various 
values of the ratio m / Mq . The ratio of linear mass approximation to full result is shown in the fourth 
line. 
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We emphasize that the light quark vacuum polarization has to be fully taken into account without any 
expansion. Comparing Eq. ( ^^ with Eq. (18), we also see the equality to —1/2 of the static potential 
at zero distance, Eq. (^3|). This explicitly demonstrates the cancellation of the linear light quark mass 
terms in the total static energy (32) at order for the case where the expansion in m is allowed. In 
Table | the size of the order light quark mass corrections is displayed for the full result and the 
linear mass approximation for ?ti/Mq between 0.1 and 0.5. The ratio of the light quark mass terms 
versus the full result is also shown. For ui/Mq = 0.3 (0.5) the linear mass terms account for 116% 
(124%) of the full result. We see that the linear mass approximation is remarkably good, even for 
rather large values of m/ Mq, owing to the vr^ enhancement factor. 

As already mentioned before, the full form of the order light quark mass corrections is 
not yet known. However, we can compare the full result and the linear mass approximation for the 
corrections that can be obtained from insertions of one-loop light quark vacuum polarizations into the 
one-loop gluon line. For illustration, in Eq. ( [2^ ) the light quark mass corrections that can be obtained 
from insertions of one-loop light quark vacuum polarizations correspond to the terms that involve the 
function P. (In (J^iLive, Eq. (|39[), the order linear light quark mass corrections that can be obtained 
from one-loop light quark vacuum polarization insertions correspond to the terms that remain after 
setting ci = C2 = di = (^2 = 0.) In the relation between the heavy quark pole and the MS mass, the 
light quark mass vacuum polarization corrections at order read 
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0.1 


0.2 


0.3 


0.4 


0.5 


i-(3),vac,fun /r4 /as \3] 

"massive / IS\ ty J J 


0.689 


1.521 


2.419 


3.352 


4.305 


e(3),vac /li(as \3] 
"massive / L 3 V 77 ' J 


0.536 


1.072 


1.608 


2.144 


2.681 


r(3),vac / (-(3) , vac, full 
"massive / "massive 


0.78 


0.71 


0.67 


0.64 


0.62 



Table 2: Comparison of the full result (second line) and the linear mass approximation (third line) 
for the order light quark mass corrections in the heavy quark pole-MS mass relation for fi = m, 
for various values of the ratio m / Mq . The ratio of linear mass approximation to full result is shown 
in the fourth line. 
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and the function P is defined in Eq. (|T5|). In the linear mass approximation, the functions fp , fj 
and fpp read (r = m/ Mq) 



37r2 



linear in m 



f(i) 



(r) 



-ln(r)+3 1n(2)-- 



vr^ r 



Fully analytic results for fp'^^ fp^'' and fpp can be found in App. In Table |2| the size of the order 
light quark vacuum polarization corrections is displayed for the full result and for the linear mass 
approximation, for m/ Mq between 0.1 and 0.5 and fj, = m. The ratio of the linear approximation to 
the full result is also shown. From the first line of Eq. (^), it is clear why one expects the linear mass 
approximation to become better at higher orders: higher orders of perturbation theory correspond to 
either larger powers of the subtracted vacuum polarization function P, or to larger powers of Ing^. 
Both enhance the contributions coming from the low momentum region g ~ m with respect to the 
contributions from q ~ Mq. We note that higher powers of P are not suppressed for small q due to the 
subtraction of the massless vacuum polarization, see Eq. (15). However, the numbers in Table || show 
that, in contrast to this argument, the linear mass approximation for the order vacuum polarization 
contributions is in fact much worse than for the order light quark mass corrections. For m/ Mq = 0.3 
(0.5) the linear mass terms only account for 67% (62%) of the complete result. Having a closer look at 
the second line of Eq. (^), however, one finds that this behaviour is caused by a cancellation between 
fp'^^ and fp^\ and the fact that the corrections to their linear mass approximations have a different 
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n 





1 


2 


3 
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5 


6 


7 


8 


9 


10 




1.17 


0.88 


1.10 


0.92 


1.07 


0.94 


1.05 


0.96 


1.04 


0.97 


1.03 


rp(n)i ,An} 
[ippjliii/ J-pp 


1.016 


1.003 


1.001 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 



Table 3: The functions fp and fpp , which describe the corrections coming from hght quark vacuum 
polarization insertions into the one-loop gluon line of the heavy quark pole-MS mass relation, for 
various values of n and for m/Mq = 0.3. 



^/[GeV] 


1.0 


1.5 


2.0 


2.5 


3.0 


3.5 


4.0 


4.5 


5.0 


<-(2),full e(3), vac, full 
"massive ^massive 


0.0750 


0.0488 


0.0383 


0.0324 


0.0286 


0.0259 


0.0239 


0.0223 


0.0210 


^(2) 1 r(3),vac 
'-^massive ' '-^massive 


0.0607 


0.0455 


0.0375 


0.0325 


0.0292 


0.0268 


0.0249 


0.0234 


0.0222 


linear app./full 


0.81 


0.93 


0.98 


1.00 


1.02 


1.03 


1.04 


1.05 


1.05 



Table 4: Comparison of the full result (second line) and the linear mass approximation (third line) for 
the sum of the order and light quark mass corrections in the bottom quark pole-MS mass relation 
for m = 1.5 GeV, Mq = 4.2 GeV and ar^^z) = 0.118, for various values of the renormalization 
scale fi. The ratio of linear mass approximation to full result is shown in the fourth line. 



sign. Individually, on the other hand, the linear mass approximations are indeed improved for larger 
powers of the logarithm and of P. This is demonstrated in Table |^, where the ratios of linear mass 
approximation to full result for fp""^ and fpp^ are shown for m/ Mq = 0.3 and n = 0, . . . , 10. 

Despite the rather good approximation provided by the linear mass terms at large orders, 
the fact remains that they only account for about 70% of the order light quark mass vacuum 
polarization corrections, for values of m/Mq that are relevant to the charm mass effects in the bottom 
pole-MS mass relation. At this point it is instructive to compare the sum of the order a1 light quark 
mass corrections and the order light quark mass vacuum polarization corrections displayed in 
Eqs. (^) and (^), (^missrv' + '^mis'swc'"") with the corresponding linear mass approximation of the sum, 
<^malsivc + '^malivc • Tablc |^ shows thc full result (second line) , the linear mass approximation (third line) 
and the ratio (fourth line) for m = 1.5 GeV, = 4.2, af'^^\Mz) = 0.118 and n = 1-5 GeV, where 
^ is the scale used in the strong coupling. For the strong coupling we used four- loop running |36] 
and three-loop matching conditions ||3^ at the five- four quark flavour threshold at /i = M^. It is an 
interesting observation that, although the linear approximation differs from the full result by +18% 
at order a1 and —35% at order a\ (at n = m), the difference amounts to only —7% for the sum. 
Therefore, instead of using the full result at order a1 and the linear mass approximation at order a^, 
it seems more advantageous to employ the linear mass approximation at order and a^. [We remind 
the reader that the order terms are as large as the order ones, see Sec. |2.3| .] This is the reason 
why, in this work, we will employ only the linear approximation for the charm quark mass corrections 
in the bottom pole-MS mass relation, as displayed in Eq. (p^). We emphasize that this decision is 
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only based on the examination of the order Ught quark mass corrections coming from vacuum 
polarization insertions. They account for about 40% (// = m) to 70% (// = 3m) of the full result in 
the linear mass approximation. A calculation of the full order light quark mass corrections would 
be useful to remove this source of uncertainty for the charm mass corrections determined in this work. 

4 Light Quark Mass Corrections in the 
Heavy Quark Pole— IS Mass Relation 

The heavy quark pole mass is quite a convenient parameter to carry out bound state calculations 
because, in the pole mass scheme, the Schrodinger equation, which describes the non-relativistic 
dynamics of a perturbative heavy-quark-antiquark pair at NNLO, has the simple form shown in 
Eq. (2). However, as has been noted before, the heavy quark pole mass parameter should not be 
employed in an actual analysis for the numerical determination of the heavy quark mass. In practice, 
in the pole mass scheme, the determination of the heavy quark mass parameter is spoiled by a large 
correlation to the value of the strong coupling, large perturbative corrections and large renormalization 
scale variations. Thus it is advantageous to keep the pole mass parameters for the actual bound state 
calculations, but to eliminate it later in favour of a better defined mass with a reduced sensitivity to 
small momenta and without an ambiguity of order Aqcd- Such mass definitions are called "short- 
distance masses". The well-known MS mass definition is the prototype of a short-distance mass. 
However, for the description of the heavy-quark-antiquark bound state dynamics through Eq. (|l]), the 
MS mass definition turns out to be a rather bad choice because it breaks the non-relativistic power 
counting. This breakdown can be visualized in Eq. (|^): after eliminating the heavy quark pole mass 
in favour of the MS mass 

M-'^ Mq(Mq) + ^(^)Mq(Mq) + ... , (52) 

the additional correction terms coming from the external energy term E = y/s — 2Mq°^" dominate all 
other terms, which are proportional to a1 and higher powers of ag- In practice, if we extract the heavy 
quark MS mass with this method, while taking care of the proper cancellation of the contributions that 
are linearly sensitive to small momenta, the results have a large correlation with the choice of the strong 
coupling and a large renormalization scale dependence, which is comparable to the pole mass case. 
(See Ref. |^8| for a comparison of the heavy quark pole, MS and IS mass definition in the framework 
of bottom and top mass determinations from bound state calculations.) From the conceptual point 
of view the heavy quark MS is not suited to the quark-antiquark bound state problem, because it is 
a short-distance mass genuinely designed for processes where the heavy quark is very off-shell. This 
means that it is a good choice for high energy processes (such as heavy quark production at high 
energies), or for processes where the external energies are much smaller than the heavy quark mass 
(such as virtual heavy quark corrections in low energy processes). A short-distance mass definition 
that is suited to processes where the heavy quark is very close to its mass-shell (i.e. for situations 
where - <^ M^) is the IS mass. The IS mass was introduced in Refs. |^, ^] to address the 
previously mentioned issues for bottom and top quark mass extractions from the T sum rules and from 
top-antitop quark pair production close to threshold at a future Linear Collider. Other heavy quark 
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mass definitions with similar properties were proposed in Ref. ||4^ (kinetic mass) and in Ref. (PS 
mass) . 

The heavy quark IS mass is defined as half of the perturbative contribution to the mass of a 
J^*-" = 1 , '^Si quark-antiquark bound state, assuming that the heavy quark is stable. In the case 
of the bottom quark, it corresponds to half of the perturbative contributions of the mass of the T(15) 
meson. We emphasize that the IS mass is a priori a purely formal parameter. Its relation to the heavy 
quark pole mass, and any other mass definition that is defined perturbatively, is determined using the 
non-relativistic power counting for perturbative heavy-quark-antiquark systems (i.e. assuming that 
the scale hierarchy Mq » MqV ^ Mqv'^ » Aqcd is realized) regardless, whether this assumption is 
valid in reality or not. For the top-antitop quark system it is certainly valid in most cases, for the 
bottom-antibottom quark system it is debatable for low radial excitations, whereas for the charm- 
anticharm system it is definitely not valid at all. Nevertheless, the IS mass could be employed in all 
cases through its perturbative relation to other mass definitions. However, the IS mass is most useful 
for systems that are truly perturbative, because for these systems it has a direct physical interpretation 
as representing the perturbative contribution of the mass of the J^*^ = 1 , ^Si heavy quarkonium 
ground state. 

In this section we present the NNLO light quark mass corrections to the relation between the 
heavy quark IS and pole masses. As in the previous section, we will assume that there are ni light 
quark species, of which one has a finite (pole) mass m, and we adopt the convention that the strong 
coupling evolves with n/ quark flavours. 



4.1 Results 

The relation between the heavy quark IS and pole masses at NNLO in the non-relativistic expansion 
including the light quark mass corrections can be parametrized as (a^ = al"'^(/x)): 
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where the subscript "massless" indicates the corrections for masslcss light quarks and the subscript 
"massive" the corrections from the light quark mass. The massless corrections are well known and 
read 
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where 

L = ln( ^ r) ■ (57) 

The one- and two- loop coefficients of the beta- function, Pq and and the constants ai and 02 are 
given in Eqs. (0)^(0) • The NNLO corrections for massless light quarks were first calculated in Ref. 
The contributions at LO, NLO and NNLO have been labelled by powers of e, and e'^, respectively, 
of the auxiliary parameter e = 1, that has already been employed for Eq. ( |34|) and that will be relevant 
in Sec. |5| when we determine the relation between the heavy quark MS and IS masses. The results at 
NLO and NNLO are obtained by starting from the known S-wave ground state solution (r = |r|), 

((>-^s{r) = (r I IS) = (lS|r) = 7r"3^i e-T^ (58) 
7 = (59) 



of the LO Schrodinger equation, 

,po,e+K"°W-2M-"'A^° )</.i5(r) = 0, (60) 



V2 



Ml 

and by using Rayleigh-Schrodinger time-independent perturbation theory to calculate the higher order 
corrections. The formal results for the light quark mass corrections at NLO and NNLO can be quickly 
written down using the Dirac notation: 
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' // J \ I c,m TTi 771 c, massless I / J V / 



where V^XZa^S^)^ '^K'^n°(^)> and (5V;™^°(r) are given in Eqs. (|), (||) and (|0D, respectively. Details 
of the calculations required by Eqs. (^ij) and (|6^) are presented in App. |A[ The final results for A|^^°;^g 
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where 



and 
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The function / is defined in Eq. (fi7\). It is instructive to consider the hmit m <C MqOs of Eqs. 
and (M), 



massive ^ \ _ / " " 
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If m « MqOs, as is the case for the charm mass effects in the bottom pole-lS mass relation, the 
expansion in the charm quark mass is completely meaningless as each term in the expansion is of 
equal size. Prom Eqs. ( [75| ) and ( [7^ ) we see that the linear light quark mass terms coincide with those 
that can be obtained from the corresponding expansion of the light quark mass corrections to the static 
potential, Eq. (^). This is expected, since the linear light quark mass terms in the static potential 
are r-independent and because all LO states form an orthogonal set. Thus, multiple insertions of the 
static potential at higher orders of Rayleigh-Schrodinger perturbation theory do not lead to any linear 
light quark mass terms. If the MS definition is used for the charm quark mass, m = m(m)[l + |(^)], 
we have to apply the replacement 

A::t°ve("i,M-'%a,) Af„tLe(^(^),M-^a,) + 6^(^)A::i,(m(m),M-^^ (77) 

in Eq. (|53|), where 

A:L1. = «^^™ = ^(£){^o(a) + l}, (78) 

and ho is defined in Eq. (p^). 

In order to implement the IS mass into the moments of the bottom-antibottom cross section 
for the sum rule analysis we need the inverse of Eq. (^3[), using the non-relativistic power counting, 
i.e. with respect to the non-relativistic expansion in as (and not in e). The result reads 



MP°'° = M^^ 1 1 + A^°(a,) + 



+ 
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A-°(a.) + AriL(M^^a.,^) + A™ (m,M^^a„^) 



4.2 Double Insertion of the NLO Static Potential 



(79) 



We already mentioned after Eq. (|76| ) that in the limit m — > the linear light quark mass terms in 
Amassi^e ouly arlsc from the single insertion of the NNLO light quark mass corrections to the static 
potential (first term on the RHS of Eq. (|62[)). This means that — if one is allowed to make the 
expansion in the light quark mass — the light quark mass corrections coming from multiple insertions 
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are subleading with respect to those arising from single insertions. For the charm mass corrections, 
however, an expansion in the charm mass is a priori not possible because (1/t) ~ My^Og ~ "Ticharm- 
Thus one can not necessarily argue that charm mass corrections arising from multiple insertions are 
smaller than single insertion ones, at least in low orders of perturbation theory. In high orders of 
perturbation theory, on the other hand, charm mass corrections coming from multiple insertions are 
subleading because, as already mentioned, the contributions from small momenta enhance the linear 
mass terms with respect to higher powers of the charm mass. 

In this section we examine the size of the NNLO charm mass corrections in the bottom IS 
mass, AJ^^j;;°^, coming from double insertions of the NLO static potential (second and third terms on 
the RHS of Eq. (|62[)). This will be important for the calculation of the NNLO charm mass corrections 
to the sum rules in Sec. ^, where we neglect the double insertion contributions. From the results given 
in App. the NNLO light quark mass corrections coming from double insertions read 




In Fig. I we have displayed the ratio AT.ss?.c,J ^TJ£le for M^°"= = 4.8 GeV, af^(Mz) = 0.118 and 
1 GeV< n < 5 GeV for m = 0.1 GeV (dotted line), 0.5 GeV (dash-dotted hue), 1.0 GeV (dashed 
line), 1.5 GeV (long-dash-dotted line) and 2.0 GeV (solid line). For the charm case (m ~ 1.5 GeV) 
we see that the double-insertion contributions amount to about 1% for ~ 1.5 GeV and increase to 
about 10% for ~ 5 GeV. This shows that the NNLO corrections already are considerably influenced 
by momenta smaller than (1/t) ~ M^a^j and that the single- insertion contributions give the dominant 
contributions to the NNLO light quark mass corrections even for the charm quark case. Thus, in 
view of the fact that the use of the linear mass approximation for the charm mass corrections in the 
bottom pole-MS mass relation introduces a relative uncertainty of around 10% in the charm mass 
corrections, it would be consistent to neglect all double insertions at the NNLO level. This is useful 
because the calculation of the double insertions is considerably more time-consuming than that of the 
single insertions. For this reason we believe that it is save, at the level of 10% accuracy for the charm 
mass corrections, to neglect the double-insertion contributions for the NNLO charm mass corrections 
to the moments of the T sum rules. 
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Figure 2: The ratio of the double-insertion contributions to the full result for the NNLO light quark 
mass corrections in the bottom quark IS-pole mass relation, A™';;°^^/A™';;°^ for Mq°'° = 4.8 GeV, 
af\Adz) = 0.118 and m = 0.1 GeV (dotted line), 0.5 GeV (dash-dotted line), 1.0 GeV (dashed line), 
1.5 GeV (long-dash-dotted line) and 2.0 GeV (solid line) plotted over the renormalization scale ^i. 



5 Heavy Quark MS- IS Mass Relation and Upsilon Expansion 

In this work we use the bottom quark IS mass as the mass definition that is extracted directly from the 
T sum rule analysis. The bottom quark MS mass is then determined in a second step, using Eqs. ( |3^ 
and (]5^). However, some care has to be taken in the determination of the perturbative relation 
between the heavy quark MS and IS masses to ensure the proper cancellation of the large corrections 



that are associated with the linear sensitivity to small momenta of the pole mass in Eqs. (34) and 
(^). Care is needed because Eq. (|53|), which gives the IS mass in terms of the pole mass, represents 
an expansion in the framework of the non-relativistic power counting, where as ~ v. Equation (|3^), 
on the other hand, which gives the pole mass in terms of the MS mass, represents a usual Feynman 
diagram expansion in the number of loops, where the power of Os corresponds to the number of strong 
coupling constants that occur in the corresponding diagrams. This means that Eqs. (|3^ and (|53|) 
cannot simply be combined by using the expansion in terms of ag. Rather, a modified perturbative 
expansion, the "upsilon expansion" ||2^, has to be employed. The upsilon expansion gives the following 
prescription: the LO, NLO and NNLO contributions in Eq. (^) are of order e, and e^, respectively, 
of the auxiliary parameter e = 1. In Eq. (p^), on the other hand, the terms of order a" are of order 
e" in the upsilon expansion. When the MS mass is then expressed in terms of the IS mass one, must 
use the expansion in e. This means that, in each order of e, different orders in are mixed. This also 
means that a determination of the heavy quark IS mass at NNLO in the non-relativistic expansion 
corresponds to an order determination of the heavy quark MS mass. If the mechanism incorporated 
in the upsilon expansion is not used, the large corrections mentioned previously can survive and lead 
to systematic errors. The upsilon expansion prescription can be understood from the asymptotic large 
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order behaviour of the "massless" N"LO perturbative coefficients c„ in Eq. (p_ 



/x(n + l)!(2/3o 



(81) 



n>l 



that arises from the contributions that are hnearly sensitive to small momenta and that dominate c„ 
for large n. This relation comes from the fact that terms that are linearly sensitive to small momenta in 
c„ involve powers of the logarithmic term L = \b.[^ / C pasMo) , which exponentiate at higher orders, 
XliLo-^V^' ~ exp(L) = ^/CpdsMQ, and effectively cancel one power of (2^. Another way to 
visualize the mechanism behind the upsilon expansion is provided by the light quark mass corrections 
themselves, as they provide a natural probe for the linear sensitivity to small momenta: the linear light 
quark mass terms in A|^^°j^^ and AJ^^j;;°^ in Eqs. (75) and ( [7^ ) are of order and a^, respectively; 
^mas°ivc ^mas'l^°e! ou thc othcr haud, are of order and af, respectively, in the non-relativistic 
power counting. The origin of this mismatch is that the non-relativistic dynamics of a perturbative 
heavy-quark-antiquark system contains scales with powers of Og such as the inverse Bohr radius MqUs 
as dynamical scales. These scales have to be accounted for by the non-relativistic power counting. 
The question of sensitivity to small momenta and the corresponding counting of orders, on the other 
hand, is only relevant to scales that are much smaller than any dynamical one, and therefore not tied 
to the non-relativistic power counting. 



5.1 Results 



Combining Eqs. (34) and (53) using the upsilon expansion up to order e we arrive at the following 



result for the relation between the heavy quark MS and the IS masses (d^ = a^r^\M^)), 
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(82) 



where we have separated the corrections for massless light quarks (subscript "massless") and those 
coming from the light quark mass (subscript "massive"). Similar to what we did in the previous 
sections, we consider n; light quarks, of which one has a MS mass m(m). The individual results for 
the corrections for massless light quarks read 

<5«, 
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(83) 
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(85) 



^ We emphasize that the issue of infrared sensitivity to small momenta, usually also referred to as the "infrared 
renormalon problem", is a purely perturbative one. The typical hadronic scale Aqcd does not arise as a relevant 
dynamical scale, but only as a dimensionful parameter that parametrizes the renormalon ambiguities in the perturbative 
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where 









^NLO 






^NNLO 


masslcss \ 




5(1) 








= (5^*^ (a ) 

— '^masslcss V'-l'Sy ) 


i = 2,3 



(86) 
(87) 



(90) 

The formulae for the functions on the RHS of Eqs. (^)-(|90|) can be found in Eqs. (|35|) -(|37|) and 
(^)~(|56|). The individual results for the light quark mass corrections read 
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The formulae for the functions on the RHS of Eqs. (|93|)- (l97|) can be found in Eqs. (|3S|), (|39|), (|63|), 
(|6^) and (^8|) . The results for massless light quarks have already been explicitly presented in Ref. |]38|| . 
It is worth mentioning that, in the limit m(m) 0, A^f^ in Eq. ( |92[ ) still contains a small linear light 
quark mass term, m(rn)af /27. This term arises because the heavy quark IS mass receives, through 
the kinetic energy term — V^/Mq"'" in the Schrodinger equation (|l|), an additional dependence on the 
pole mass. This does not affect the cancellation of linear light quark mass terms in the total static 
energy. The associated ambiguity exists, regardless of which of the currently known short-distance 
mass definitions is employed in an extraction of the bottom quark mass parameter from non-relativistic 
bottom-antibottom observables. However, the ambiguity is suppressed by A^° = CpO^/S = 0.02 
with respect to the ambiguity of the pole mass definition, i.e. it is smaller than the order Aqqy)/Mq 
ambiguity for the bottom quark case. It can therefore be ignored for all practical purposes. 

For simplicity we have displayed Eq. (|82|) for the strong coupling at the scale ^ = M^. 
Equation (|82|) can be generalized to arbitrary scales with the relation (a 
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Figure 3: Charm mass corrections in the relation between the bottom MS mass Mb(Mb) and the 
bottom IS mass. Figure (a) displays the scale dependence for M^^ = 4.7 GeV and ai^^(M^) = 0.118 
and using m(m) = 0.1 GeV (dotted line), 1.0 GeV (dash-dotted line), 1.5 GeV (dashed line) and 
2.0 GeV (solid line) for the MS charm quark mass. Figure (b) displays the dependence on m(m) for 
fj, = 1.5 GeV (dotted line), 3 GeV (dashed line) and 5 GeV (solid line). The other parameters a chosen 
as in Fig. (a). 



We note that, when we examine the scale dependence of the light quark mass corrections in the 
following subsection, we leave the scales of the quark masses unchanged. We also note that from 
now on we will designate to the terms of order e" in Eq. (|8^ ) by the expression "order a^". This 
accounts for the fact that a determination of the heavy quark IS mass at NNLO in the non-relativistic 
expansion corresponds to an order determination of the heavy quark MS mass, as far as the proper 
cancellation of the linear infrared-sensitive contributions is concerned. 



5.2 Brief Examination 



It is instructive to examine the size of the light quark mass corrections in Eq. ( p2[ ) for the case of the 
bottom quark (n; = 4) taking into account the mass of the charm quark and treating the other light 
quarks as massless. For M^^ 



4.7 GeV, m{m) = 1.5 GeV and ai^V = 4.7 GeV) = 0.216 we obtain 
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0.030 + 0.0049„ 



(99) 



The corrections coming from the finite charm quark mass are indicated by the subscript "m". The 
shift in Mb (Mb) caused by the charm mass corrections is —7.2 MeV at order (e^) and —4.9 MeV 
at order (e^). If we neglect the NNLO double-insertion corrections in A^^j;^°^ in Eq. (85), we 
obtain —6.1 MeV at order a^. It is quite instructive to confront the perturbative series of the bottom 
pole-MS and the bottom IS-pole mass relations in Eqs. ( [3^ ) and (^3|) with the result for the MS~1S 
mass relation in Eq. (|9|). For M^"^" = 4.9 GeV and Mb (Mb) = 4.2 GeV and the same values for a^, 
/u and m(m) as for Eq. (pW), we obtain: 
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Table 5: The scale dependence of the "massless" (no subscript), "massive" (subscript "m") and 
the large-/3o (subscript "/So") corrections in the bottom MS-IS mass relation for m(m) = 1.5 GeV, 
M^^ = 4.7 GeV and a^^\Mz) = 0.118. The superscript "(i)" corresponds to corrections at order a* 
(or e' in the upsilon expansion). 



Clearly, the bottom MS-IS mass relation has a much better convergence for the corrections for massless 
light quarks as well as for the charm mass corrections. As expected, we find that the cancellation of 
the linear charm mass term is more efficient at order (e^) than at order (e^). In Fig. ^a the 
scale dependence of the charm mass corrections M^^(A^) + Af^^) is displayed for M^^ = 4.7 GeV and 
af\Mz) = 0.118 for m(m) = 0.1 GeV (dotted line), 1.0 GeV (dash-dotted line), 1.5 GeV (dashed 
line) and 2.0 GeV (solid line). For the strong coupling we have employed four-loop running and three- 
loop matching conditions at the five-four quark flavour threshold. In Table |5| the scale dependence of 
the "massless" and "massive" corrections in Eq. (p^) are shown for rn{m) = 1.5 GeV, M^^ = 4.7 GeV 
and Us (Mz) = 0.118, separately for each order. We see that the convergence of the charm mass 
corrections is best for ~ m(m), whereas for the massless corrections this happens for fi ~ 2.5 GeV. 
This reflects that the characteristic scales for the two types of corrections are slightly different. This 
is not unexpected as the relevant dynamical scales for the charm mass corrections are the charm mass 
and the inverse Bohr radius My^Os', the massless corrections, on the other hand, involve and My^ag 
as relevant dynamical scales. The convergence of the charm mass corrections shown in Table |5| clearly 
shows that they are well under control. For our final determination of the bottom quark MS mass in 
Sec. 0, we will use renormalization scales between 1.5 and 7 GeV. For the determination of the central 
value of the bottom MS mass we will use ^ = M^^. 

In Fig. ^b, finally, the size of the charm mass correction M^^(A^) -|- Aj^^) is plotted over m(m) 
for = 4.7 GeV, al^^(Mz) = 0.118 and n = 1.5 GeV (dotted hue), 3 GeV (dashed line) and 5 GeV 
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(solid line) to illustrate the impact of other light quark masses on Mb(Mb). Figures ^ and b show 
that for m(m) = 0.1 GeV the corrections amount to much less than 1 MeV. Thus, the effects of the 
masses from light quarks other than the charm can be safely neglected. 



5.3 Estimate of the Order Charm Mass Effects 

The results shown in the 9th and the 10th row of Table |5| reveal a quite interesting property of 
the ratio between the full corrections and the massless ones at order and a^. For ^ around 4 
to 5 GeV the ratios at order and are approximately equal to (/3o"'~^V/'o"'~^^) ~ ^-^^ 
(/?q"'~^^//3q"''~^'')^ = 1.17, respectively. We already mentioned that the inclusion of the charm mass 
leads to an effective decoupling of the charm quark corrections at large orders. Thus, because the 
order a""*"^ term in the MS-IS mass relation in Eq. (]8^) contains (]q, the value {P^'~^^ / 
for the ratio might not be unexpected as the asymptotic behaviour for large orders of perturbation 
theory. However, we emphasize that if ones employs ai"'~^^ rather than ai"'''~^\ as we do in this 
work, the decoupling is not explicit. (See the results for the light quark mass corrections for the static 
potential discussed in Sec. |2|.) In addition, because the series in the bottom MS-IS mass relation 
in Eq. ( |82| ) involves the cancellation of the leading infrared-sensitive contributions, also subleading 
contributions, which involve for example lower powers of /5o or higher order contributions in the 
non-relativistic expansion, will contribute to the dominant asymptotic behaviour of the perturbative 
coefficients in Eq. (^%. Thus the observation above is certainly a coincidence, particularly the one at 



order a^. Nevertheless, the ratio (/3q"''~^^//5q"' "'''')" might be used as a rough estimate for the large 
order behaviour of the ratio of full to massless corrections in the MS-IS mass relation. It should 
therefore be possible to estimate the size of the order charm mass corrections from the order 
corrections for massless light quarks, A^^^. At the present stage the full result for A^^^ is not yet 
known. However, it can be rather easily calculated in the large-/3o approximation because for the N^LO 
corrections to the pole-lS mass relation this only amounts to the insertion of one-loop massless quark 
vacuum polarizations into the LO static potential. The order large-/9o corrections to the pole-MS 
mass relation are already known from Ref. |41|, whereas the N^LO corrections to the pole~lS mass 
relation in the large-/3o approximation are calculated in App. M. The formula for A^^^ in the large-/3o 



approximation is derived in App. |C| and reads ( 
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(102) 



only involves corrections proportional to /3o 
(i = 1,2,3,4) for m(m) = 1.5 GeV, M^^^ 



where L = ln(l/Ciras). Equation ( |102| ) can be generalized to arbitrary renormalization scales in 
the strong coupling using the one-loop running of the strong coupling, in analogy to Eq. (p^), which 

We have displayed the scale dependence of M^^^A^l| 
4.7 GeV and af\Mz) = 0.118 in the 11th to the 
14th rows of Table |5|. The values for very small choices of /i are meaningless. As indicated in the 
discussion about the ratio (/Jo^V/^o"^^)" the beginning of this subsection, we see that the large- /Jq 
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approximation does certainly not provide a very good estimate of the true corrections. At order 
and the large-/3o result overestimates the true correction. For fi ^ 3 GeV the difference is between 
50 and 100%. This demonstrates that the large- /3o approximation does in general not give more than 
an order of magnitude estimate if it is applied to a perturbative expansion that does not contain 
a linear sensitivity to small momenta, as is the case in the MS-IS mass relation. (In the pole-MS 
mass relation, which contains a linear sensitivity to small momenta, the large-/3o approximation is 
considerably better, see e.g. Ref. On the other hand, the large-/?o approximation seems to 

indicate quite a good convergence for larger values of Although we cannot prove that this property 
allows for any conclusion about the size of the true expression for Ml^A^^\ we believe that this is 
likely to be the case^. If the relative size of the large- /3o corrections roughly reflects the relative size of 
the true corrections, one can expect that M^^A^'^^ amounts to around —10 to —15 MeV for around 
5 GeV. Using the value (Po^/f^Q^^)^ = 1.26 as an estimate for the ratio (A^ A(f))/A('^), this results 
in A^) = —3 to —4 MeV for fi around 5 GeV. This estimate seems to be consistent with the sum 
M^^(A[^) -|- A^';')) at the scale /i = 1.5 GeV, where we find the best convergence of the charm mass 
corrections. 

To conclude the somewhat speculative discussion of this subsection, we note that the main 
outcome is that the large- /?o approximation does not provide a particularly good estimate of higher 
order corrections in the MS-IS mass relation. This is because subleading infrared-sensitive contribu- 
tions would have to be taken into account to determine the dominant asymptotic behaviour of the 
perturbative relation of short-distance masses. As far as our estimate for the order charm mass 
corrections in the bottom MS-IS mass relation is concerned, we believe that it is more likely to be 
correct than the estimate for the order massless corrections, as the light quark mass corrections at 
order do not yet contain any ultrasoft contributions. 



6 Light Quark Mass Effects in the T Sum Rules 

In this section we determine the corrections coming from the finite charm quark mass to the large-n 
moments of the bottom-antibottom quark cross section in e~^e~ annihilation at NNLO in the non- 



relativistic expansion. In Sec. |6Jj we briefly review the basic concepts involved in the T sum rules. 
More details can be found in Ref. In Sec. we outline our method to calculate the moments, 
and in Sec. |6.3| we present the calculations of the charm mass corrections. Results for the calculations 
for massless light quarks are not given here; they have been given in Ref. 0. In Sec. |6^ the results 
for the charm mass corrections are examined. 

6.1 Basic Issues 

The sum rules for the T mesons start from the correlator of two electromagnetic currents of bottom 
quarks at momentum transfer q, 

U^,{q) = -i fdxe''^-^{0\Tjl{x)jt{0)\0) 



* We are aware of the fact that A''*^ contains ultrasoft "Lamb-shift" -type corrections that involve MbQ^ b^s the relevant 
dynamical scale and might lead to considerable additional contributions. Calculations of the latter have been carried 
out in Refs. However, no concrete statement can be made before all corrections in A'*' are determined. 



32 



^ lm[-i(0|rj^(<?)?''^(-<?)|0)], (103) 

where 

4(x) =b{x)j,b{x), (104) 

and the symbol b denotes the bottom quark Dirac field. The n-th moment Pn of the vacuum polar- 
ization function is defined as 



g2=0 



n!g2 \dq 

where Qh = —1/3 is the electric charge of the bottom quark. Due to causality the n-th moment P„ 
can also be written as a dispersion integral 

oo 

Pn= f ^ R{S) , (106) 



where 



CTpt S ^ 

is the total photon-mediated cross section of bottom quark-antiquark production in e+e^ annihilation 
normalized to the Born cross section for massless leptons, apt = 47ra^/3s, and s is the square of the 
centre-of-mass energy. The lower limit of the integration in Eq. ( |106| ) is set by the mass of the lowest- 
lying resonance. Assuming global duality the moments Pn can be either calculated from experimental 
data on R or theoretically using perturbative QCD. 

The experimental moments P^ are determined by using latest data on the T meson masses, 
Mks, and electronic decay widths, Fks, for /c = 1, . . . , 6. The formula for the experimental moments 
used in this work reads 

6 °° 

Pn = ^Y.^.+ I ^W(.), (108) 



(v^)bb 



and is based on the narrow width approximation for the known T resonances; ciem is the electromag- 
netic coupling at the scale 10 GeV. Because the difference in the electromagnetic coupling for the 
different T masses is negligible, we chose 10 GeV as the scale of the electromagnetic coupling for all 
resonances. The continuum cross section above the BB threshold is approximated by the constant 
Tc = 1/3, which is equal to the Born cross section for s oo, assuming a 50% uncertainty, 

rcont(s) = r,(l±0.5). (109) 

For n > 4 the continuum contribution is already suppressed sufficiently so that a more detailed 
description is not needed. For a compilation of all experimental numbers used in this work, see 
Table |. 

A reliable computation of the theoretical moments P^ based on perturbative QCD is only 



possible if the effective energy range contributing to the integration in Eq. (106) is sufficiently larger 



than Aqcd ~ 0(200-300 MeV) ||. For large values of n it can be shown that the size of this energy 
range is of order M^/n. Qualitatively this dependence on the moment parameter n can be seen from 
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nS 


Mns/[GeVl 


rns/fkeVl 


IS 


9.460 


1 32 + 04 ± 03 


2S 


10.023 


0.52 ±0.03 ±0.01 


3S 


10.355 


0.48 ±0.03 ±0.03 


4S 


10.58 


0.25 ±0.03 ±0.01 


5S 


10.87 


0.31 ±0.05 ±0.07 


6S 


11.02 


0.13 ±0.03 ±0.03 



Table 6: The experimental numbers for the T masses and electronic decay widths used for the 
calculation of the experimental moments P^. For the widths, the first error is statistical and the 



second one systematic. The errors for the partial widths of T(1S) and T(2S) are taken from Ref. |44]. 
All the other errors are estimated from the numbers presented in Ref. ||2^. For the electromagnetic 
coupling at 10 GeV and the BB threshold point we use q.~^ = a^(lOGeV) = 131.8(1 ± 0.005) and 
(V^)bb = 2 X 5.279 GeV. The small errors in the T masses and the BB threshold (y/s)^^ are neglected. 



Eq. ( 106 ): large values of n enhance the low energy contributions and suppress the high energy ones. 
This implies that n should be chosen sufficiently smaller than 15-20. This is consistent with estimates 
of non-perturbative effects in the large-n moments using the gluon condensate 46 1 : 



Dth 

1 n n 



3 

1 - (asG2)-^-^exp(-0.4CFa, 



72 



(110) 



where the P^ represent the theoretical moments obtained with pure perturbation theory. The RHS 
of Eq. ( |110| ) amounts to less than a per cent of the theoretical moments for n < 20. However, this 
result needs some careful interpretation because the estimate of non-perturbative effects based on 
an expansion in gluonic and light quark condensates is only reliable if M^,v'^ ^ Aqcd- In addition, 
the values of the dimension-6 (and higher) condensates as well as all perturbative Wilson coefficients 
for the condensates are unknown, which introduces additional systematic uncertainties. In order to 
suppress these sources of uncertainties, n should be chosen as small as possible. However, it is also 
desirable to choose n as large as possible in order to suppress the contribution from the bb continuum 
to R{s) above the BB threshold, which is rather poorly known experimentally. In other words, one 
has to choose n large enough for the bottom-antibottom quark dynamics encoded in the moments 
Pn to be non-relativistic. Because the size of the energy range contributing to the n-th moment is of 
order M^/n, the mean centre-of-mass velocity of the bottom quarks in the n-th moment is u ~ l/^/n. 
This counting rule allows for a quantitative formulation of the two requirements for the choice of 
n: theoretical reliability demands M^v^ ^ Aqcd; and dominance of the non-relativistic dynamics 
demands v 1. If both requirements are met, the bottom-antibottom quark dynamics encoded in 
the moments is perturbative, i.e. the hierarchy ^ Af^f ^ M^v"^ ^ Aqcd is valid. Only if this 
condition is satisfied, the non-relativistic power counting described after Eq. (|lj) and used throughout 
this work is truly justified. Only in this case can the non-relativistic bottom-antibottom quark pair 
be considered as a Coulombic system for which the perturbative treatment carried out in this work is 
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feasible and for which the power counting 



1 , , 

~ ~ , (111) 



used for the calculation of the moments, can be applied. In our analysis we choose 

4<n<10 (112) 

as the allowed range for the moments Pn- The upper bound was chosen to avoid the problem of 
unknown systematic non-perturbative uncertainties as much as possible. The gluon condensate con- 
tributions mentioned above amounts to less than a per mille in the moments for n < 10. This is 
an order of magnitude less than the charm mass corrections. Therefore, non-perturbative effects are 
not taken into account in this work. The lower bound, on the other hand, is a compromise born 
out of the requirement of the non-relativistic dominance and of the fact that we would like to have 
a reasonable number of moments that can be used for fitting. We also note that our estimate of the 
size of non-perturbative effects is based on the common faith that the expansion in gluon condensates 
really describes the dominant source of non-perturbative contributions. This is not necessarily the 



case, as has been pointed out for instance in Refs. |^. As we will discuss in Sec. |6^, the pertur- 
bative behaviour of the moments does not seem to comply at all the estimate of the smallness of 
non-perturbative effects based on the gluon condensate (presuming that perturbative convergence has 
anything to do with the size non-perturbative effects). We believe that it is probably the weakest 
point of the sum rule method. 

6.2 Method for the Calculation of the Moments 

In this subsection we briefly review the method used in Ref. for the calculation of the NNLO 
moments for massless light quarks assuming the pole definition for the bottom mass. This sets the 



stage for the determination of the charm mass corrections carried out in Sec. . 

At NNLO in the non-relativistic expansion, we need to calculate corrections up to order a^, 
C(s/V^ 1/n to the moments P„ with respect to the non-relativistic (LO) limit. According to the 



power counting ( |lllD we need to keep corrections of order aj, agV and v n with k — 21 = 2 in the 



non-relativistic expansion of Eq. (|106D , 



pth _ 1 f dE ( E \( E E^ \ 

where E = y/s — 2Mb°'° and -Bbind is the perturbative binding energy of the lowest-lying resonance, 
i.e. Sbind = 2{M,]^ - Mr'"). The exponential exp{-E /M^'^n) is the LO term in the non-relativistic 
expansion of ds/s"^^^ and cannot be expanded because E/M^°^''n is of order 1 according to the power 
counting of Eq. ( [1111 ). We note that the non-relativistic expansion of ds/s"'~^^ does not contain any 
NLO corrections, as these would be non-analytic functions of the energy E. Thus, because the light 
quark mass corrections start at NLO, we only need to consider the LO term in the non-relativistic 
expansion of ds/s"'^^ in the determination of the light quark mass corrections up to NNLO. The 
term -R^nlo is the bottom-antibottom quark cross section expanded at NNLO in the non-relativistic 
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expansion. It is obtained by expanding the relativistic electromagnetic bottom quark currents in 
Eq. ( 107 ) in terms of effective non-relativistic ^Si currents up to dimension 5 (i = 1,2,3): 



IKq) = ci(^Va)('z)-^(^V,(-|B)2x)(g) + ..., 



(114) 



Cl 



+ 



(115) 



where s = . The non-relativistic currents are defined in the NRQCD factorization scheme proposed 
by Lepage et al. in Refs. |48, 49 1, which separates contributions coming from momenta of order 
from non-relativistic momenta of order M^v and M^,v'^. The constants ci and C2 are short-distance 
Wilson coefficients that contain the contributions from momenta of order M^. At the Born level, 
Cl = C2 = 1. At NNLO we need ci at order a^; for C2 we do not need to calculate any corrections as 
the contributions of dimension-5 currents are already suppressed by two powers of v. From Eqs. ( |114D 
and (|115|) one obtains 



-'•'NNLO 



{E) 



Cl(/^hardi Wac) I™ Al{E , ^soft; Wac) 



{Mr'' 



c^2 



polc\4 
b ) 



3(M, 



Im 



-4.2 (-E, /isoft) 



(116) 



where Ci = c\ and 



Al = i{Q\{i,^ax) (x^aV^)lO), (117) 

A2 = \i{Q\{i,^ax){x^S{-I^Dfij) +h.c.\^). (118) 

The constant = —1/3 is the electric charge of the bottom quark. Using the NRQCD equation 
of motion for the bottom and antibottom quark fields, one can relate A2 to Ai, A2 = Ur'^EAi. 
In Eq. ( |116| ) we have indicated the dependence of the non-relativistic cross section on the renormal- 
ization scales used in Refs. @, ^J. The scale /igoft is the renormalization scale of the strong coupling 
governing the bottom-antibottom quark potential in Eq. (||) and /Xhard the renormalizations scale of 
the strong coupling in the short-distance coefficient Ci. The scale yUfac is a factorization scale that 
separates non-relativistic from hard momenta in the NRQCD factorization scheme. The scales ^)\a.rd: 
/Ltfac are irrelevant to the charm quark mass corrections determined in this work, since they contain no 
hard corrections up to NNLO. The dimension-6 current correlator ^1 is obtained at NNLO from the 
configuration space Green function of Eq. (|^) for both spatial arguments evaluated at zero distance: 

^1 = 6AfcG(0,0,S) . (119) 

(See Ref. [^] for details on the derivation of Eq. ( |119[ ).) The short-distance coefficient Ci is obtained 
by matching the NNLO cross section in Eq. ( |116| ) to the same cross section in full QCD in the (formal) 
limit CKs <C f <C 1 at the two-loop level and including terms in the velocity expansion up to NNLO. This 
method of determining the short-distance coefficient Ci at the level of the cross section (rather than the 



amplitude) is called "direct matching" |50|. The calculation of the NNLO zero-distance Green function 
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of the Schrodinger equation (|l]) proceeds in analogy to the determination of the NNLO corrections 
to the IS mass carried out in Sec. ^ via Rayleigh-Schrodinger time-independent perturbation theory. 
Parametrizing the zero-distance Green function as 



G(0,0,^) = Gc(0,0,^) +G^^°(0,0,^) +G""~'L°(0,0,S) 



(120) 



where Gc is the LO non-relativistic Coulomb Green function, the formal expressions for G^^'~' and 
^NNLo 2ero distances read (r = |r|, r' = \r'\): 

G-^°(0, 0, = - j d^r Ge(0, r, E) w(^) Gc{r, 0, E) , (121) 

G™^°(0, 0,E) = - J d'r Ge(0, r, E) Kr:Le(r-) + ^Bp(r) + ^.^(r-) + 5/7,, ) G,(r, 0, i?) 

+ |ci=^r I dVG,(0,r,S)V;-J^°,Jr)Gf(r,r^£;)y-J^°^^^^ (122) 
in configuration space representation, where SH^^i^ = — V^/4Mb°^'', 
G',{ry,E) = -3- [dnG,ir,r',E) 



Mr " k sin(7r A) _k(r+r') 
2tt ttX ^ 



dt e 



-2krt 



1 + t 



du e 



2kr'- 



1-U 



u 



(123) 



GciO,r,E) = Gc{r,0,E) 



dt e 



-2krt 



l+t 



(124) 



and 



k = \l-EMr'\ 



(125) 



A = 



Gp as M, 
2k 



pole 



(126) 



The expressions for the Coulomb Green function in Eqs. ( |123| ) and ( 124 ) have been derived in Refs. |l5l[| . 
For the NNLO contributions arising from second order Rayleigh-Schrodinger perturbation theory we 
only need the S-wave component of the Coulomb Green function, Gf , for the sum over intermediate 



states because the static potential at NLO is spin- and angle-independent. The case r < r' in Eq. (123) 
is obtained by interchanging r and r' . 

The calculation of the moments as shown in Eq. ( |113D is somewhat cumbersome, because it 
requires the separate determination of bound state {E < 0) and continuum (E > 0) contributions. In 
particular, for the bound state contributions, the binding energies and the leptonic decay rates (which 
are proportional to the modulus squared of the configuration space wave functions at the origin) of all 
resonances would need to be determined individually. A more economic way to calculate the moments 
can be achieved by deforming the path of integration in Eq. (|113| ) into the negative complex plane, as 
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Figure 4: Path of integration to calculate expression ( |127D for the theoretical moments. The dashed 
line closes the contour at infinity and does not contribute to the integration. The free constant 7 is 
chosen large enough to be safely away from the bound state poles, which are indicated by the grey dots 
on the negative energy axis. The thick grey line on the positive energy axis represents the continuum. 



shown in Fig. 0. Because the path that closes the contour at infinity does not contribute, the result 



of the integration in Eq. ( |113D can be written as 
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(127) 



where the free constant 7 is chosen much larger than the ground state binding energy E'bind- In the 
second line of Eq. ( |127D the change of variable E — > —E has been performed. We note that, in contrast 
to Eq. ( |113| ), we also need the real part of the correlators ^1 and A2 in Eq. ( |127D . The expression in 
Eq. ( |127[ ) has three advantages: (i) the bound state and continuum contributions are obtained at the 
same time, (ii) the correlators, i.e. the Green functions shown in Eq. (120), can be naively expanded 
in Us before the integration, and (iii) we can use a vast number of tables for the calculations, because 
the complex energy integration is nothing else than an inverse Laplace transform. This elegant way 
to determine the moments has been used first in Ref. 

The final result for the NNLO theoretical moments for massless light quarks can be cast into 
the form 



th 



4„+l (^Polo)2n „3/2 



c 



( ^hard A^fa 



pole 



hard. 



Atsoft A^fac 



pole ' 



M, 



pole 
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(128) 
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where 



8^3/2^3/2 I ^^r dE ( E \ E E? 



7—100 



, ^3/2 ^3/2 ^ f dE / E \ r 4 E 



(129) 



7~ioo 

The expUcit results for Qn^i and ^>n,2 can be found in Ref. 0. We note that Qn,i and Qn,2 are di- 
mensionless and that Qn^i depends only logarithmically on the bottom pole mass at NLO and NNLO 
through the dependence of the strong coupling on the inverse Bohr radius. We also note that only 
Qn,i receives light quark mass corrections, since gn,2 contains only NNLO contributions. 



6.3 Light Quark Mass Corrections to the Moments 

The light quark mass corrections to the moments are determined following the lines of the calculation 
for massless light quarks described in the previous subsection. The insertions of the NLO and NNLO 
static potentials that are needed are in complete analogy to the calculation of the heavy quark IS 
mass presented in Sec. ^ 



Light Quark Mass Corrections at NLO 

The NLO light quark mass corrections to the zero-distance Green function read 



G:^''{0,0,E) = - J d^rGMr,E)SV:^^^°ir)G,ir,0,E)^ 



(131) 



where the formulae for SV^^° and are given in Eqs. (|l8| ) and (124), respectively; SV^^l^'^{r) contains a 
contribution from the massive light quark vacuum polarization function, called "massive" corrections in 
the following, and a contribution from the subtraction of the massive light quark vacuum polarization 
function in the limit m — > 0, called "massless" correction in the following. This subtraction is a 
consequence of our convention that rii quark species contribute to the evolution of the strong coupling. 
For the massless corrections we can recycle the results for massless light quarks given in Ref. 0. The 
massive corrections to the zero-distance Green function involve the term (a^ = ai^\/isoft)): 



J d^rG,{0,r,E) 



Gc{r,0,E) 



(ix/(x)e-2"^™ 

_o 00 00 

C.a. ^-^^^ Jd.m jdt JduJdrre-nH^^^^^)^m.,^ (i±^)' 
1 000 



t u 



[l + t + uf{k + a) 



n\2 ' 



(132) 
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where 



m X 



(133) 
(134) 



and k and A are defined in Eqs. (|125| ) and ( |126| ), respectively. In the last line of Eq. (132) we have 
expanded in owing to the integration path in the negative complex energy plane shown in Eq. ( |127D . 
To determine the contributions of Eq. ( |132| ) to the light quark mass corrections of the moments Qyi \ , 
we use the following inverse Laplace transform [^]: 



7+ioo 



27rz 



Z^-P 2 



7— loo 



{z + Af 



n2 gQ{p,A n) 



where 



1 + 2X 2^{l + X) 



r(i±£,x) Tilx) 
r(i±£) r(|) 



The final result for the NLO light quark mass corrections to the moments reads 



c m,NLO 
(>enA 



1 p=i 



p=i 
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m 



Mr" ' 
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00 00 00 



dt du dx f{x) 







p\ (1 + t + 'U)2 
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(l+t)(l+n) 
t u 



(135) 



(136) 



(137) 

(138) 
(139) 



9 9 



(140) 



The first term in the curly brackets on the RHS of Eq. ( |137D comes from the massive corrections 
and the second from the massless ones. The constants Wp and Wp have already been calculated in 
Ref. their expressions can be found in App. ^. The function go has been calculated numerically 
using standard variable transforms to treat the logarithmic singularities. We note that the numerical 
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determination of gQ represents a considerable effort. For given values of p and n we have calculated qq 
for r = 0.01, 0.2, 0.3, 0.4 and 0.5, and constructed an interpolation function to determine the values 
for arbitrary values of r. This was done in order to speed up the computations that arise during the 
fitting procedure. We checked that the relative deviation of the interpolation from the exact result is 
at the per cent level. We also note that the terms of the sum over p only forms a convergent series 
for p larger than 10. The larger the value of n, the more terms in p have to be taken into account. 
For n < 20 it is sufficient to carry out the sum up to p = 20, which is what we have done is our 
numerical evaluation of the moments. The issues just mentioned also apply to the NNLO light quark 
mass corrections to gn,i- 



Light Quark Mass Corrections at NNLO — Single Insertions 

The NNLO light quark mass corrections to the zero-distance Green function read 

<5Gr^°(0,0,i?) = - I d3rG,(0,r,i?)5KT°(r)G,(r,0,i?) 

+ Jd^r j d^r'G,{<d,r,E)5V:^^°{r) G!{r,r' , E) 6V::^°{r') G,{r' ,0, E) 

+ ^Jd'r J d^r' GM r, E) 5V:^^^{r) G^r, r' , E) , Jr') G,{r\ 0, E) . (141) 

Let us first consider the contributions from the single insertion of the NNLO light quark mass cor- 
rections to the static potential at first order Rayleigh-Schrodinger perturbation theory. As for the 
NLO calculation, we distinguish between "massive" and "massless" corrections and only give details 



for the massive contributions. Apart from expression (132) there is one more term coming from the 
light quark mass corrections to the NNLO static potential that is relevant for the massive corrections: 
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To determine the corrections to the moments gn,i we need the following inverse Laplace transforms |]& 
in addition to Eq. ( [1351 ): 
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where 



gi{p,X) = -2 



g2ip,X) = -2 



r(i; 



+ 1 



(146) 



+ e^X^-2 (3-P + 2X) 1 



(147) 



2 / ^ V 2 

The final result for the single insertion NNLO light quark mass corrections to the moments reads 
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The constants w^, and Wp have aheady been calculated in Ref. Q; their expressions can be 



found in App. ^. The three terms in the curly brackets on the RHS of Eq. (148) originate from the 
corresponding three terms given in Eq. (|30|), which all vanish individually for m ^ 0. If the MS 
definition is used for the charm quark mass, m = m(m)[l + |(^)], we get an additional contribution 
to the moments at NNLO through the replacement 
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where 
and 
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Light Quark Mass Corrections at NNLO — Double Insertions 



The second and third terms on the RHS of Eq. (141) come from double insertions of the NLO static 
potential. The calculation of these contributions is much more involved than that of the ones coming 
from single insertions and will not be presented in this work. On the other hand, it is not necessary to 
include them in the present analysis. We have shown that the NNLO double-insertion contributions 
in the heavy quark pole-lS mass relation are suppressed with respect to the single-insertion ones; for 
the NNLO charm mass corrections, the double-insertion contributions amount to about 10% of the 



full NNLO result. As was pointed out in Sec. 4^, this happens because at NNLO the corrections are 
dominated by momenta low enough for the linear (and r-independent) charm quark mass terms to 
dominate in the static potential. The same is expected for the double-insertion contributions of the 
NNLO charm quark mass corrections to the moments. Thus, the double-insertion contributions are of 
the same order as the terms oc {m^^^^^/M^)^ , n > 2, in the bottom quark pole-MS mass relation, which 
are already neglected in this work. For the reasons mentioned here, we neglect the double-insertion 
contributions. 



Implementation of the IS Mass 



The moments displayed in Eq. ( |12^ ), including the charm quark mass corrections to Qn,i in Eqs. (137) 
and ( |14^ ), still depend on the bottom quark pole mass through the global factor (Mr'°)"^" and 
through logarithms of ratios with the renormalization scales. The ratios of the pole mass with the 
renormalization scales originate from the running of the strong coupling and from anomalous dimen- 
sions of the ^Si NRQCD bottom quark currents. These ratios only arise at NLO and NNLO in Qn^i 
and Ci. They do not lead to any modifications in the IS mass scheme, because the difference between 
the bottom IS and the pole masses is of order ~ ~ 1/n in the non-relativistic power counting. 
The only relevant modification that arises from the implementation of the bottom IS mass comes 
from the global factor (M^"'")"^". Using Eq. (|79|) and the scaling ~ 1/n, the global factor has the 



following form in the IS mass scheme at NNLO in the non-relativistic expansion (a^ = ctf^ {iXsoit))'- 
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- 2 ( AriL(K^ a., ^soft) + Ari°.(m(m), M^, a„ ^soft) 

-2<5«(a,)A:'°(m(m),M,^^a,)] |, (160) 

where the A's have been given in Eqs. (|35|), (|5^-(^), (|63|), (|^) and (|7^). The second and third terms 
in the curly brackets on the RHS of Eq. ( |16[1| ) are the NLO and NNLO contributions, respectively. 
We emphasize that using Eq. ( |160D for the global factor l/(Mb°''')^"' in Eq. ( |128|) is equivalent to 
implementing the bottom IS mass via Eq. ( [79| ) into the Schrodinger equation (||) before any calculation 
is carried out. The renormalization scale on the RHS of Eq. ( |16[)| ) is /Ugoft, which governs the strong 
coupling in the static potential. This is because the corrections on the RHS of Eq. ( |160| ) only involve 
non-relativistic momenta. From the technical point of view, they are supposed to compensate for the 
bad large order behaviour of the static potential. For the same reason we have to consistently expand 
the LO, NLO and NNLO corrections in Eq. (|16C| ) with the LO, NLO and NNLO corrections in Qn,i- 
Aside from the issue of the cancellation of the large high order corrections that are associated with 
the bottom quark pole mass definition, it is easy to understand how the IS mass definition reduces 
the correlation of the moments (and also of the mass determination) to the strong coupling and 
the dependence on ^soft^ in the IS scheme the exponential energy-enhancement of the bottom- 
antibottom ground state that dominates the moments for large values of n, is eliminated. A detailed 
discussion of this effect has been given in Ref. Q. 

Light Quark Mass Corrections to the Short-Distance Coefficient Ci 

The short-distance coefhcient Ci contains the effects of the bottom-antibottom quark production 
process that arise from momenta of order M^. This has two consequences: the first is that the light 
quark mass corrections can be naively expanded in m? without the emergence of any terms that are 
non-analytic in m?. Thus, Ci does not contain any linear dependence on a light quark mass. The 
second consequence is that one can expand in the light quark mass because m/M^ is small. Therefore, 
the dominant light quark mass effect in Ci is proportional to (as/7r)^(m/Mb)^, where the scale in 
the strong coupling is of order M^. For the charm quark this correction is at the per mille level and 
subleading. The charm mass corrections to Ci are of the same order as the m? /M^ corrections in 
the pole-MS mass relation that have been neglected in this work. The light quark mass corrections 
to Ci are therefore also neglected in this work. We note that the absence of any linear light quark 



mass corrections in Ci has been proved by explicit calculation in Ref. [d3|. It was shown, for any 
energy at order a^, that the total cross section R, Eq. ( |107] ), does not contain any linear light quark 
mass corrections if a short-distance definition is employed for the bottom quark mass. Thus, because 
Ci = 1 (and is mass independent) at the Born approximation, it cannot contain linear light quark 
mass corrections at order ai. 



6.4 Brief Examination 

Before turning to the statistical analysis it is useful to examine the size and behaviour of the light 
quark mass corrections to the moments for various choices of the light quark mass. This will be 
helpful for the interpretation of the result that is obtained for the bottom IS mass from the sum rule 
analysis carried out in the next section. Detailed discussions on the behaviour of the moments in the 
pole and the IS mass schemes for the case of massless light quarks can be found in Refs. 01. The 
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0.698 
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Table 7: The theoretical moments P*^ at LO, NLO and NNLO for M^"" = 4.7 GeV, ay' {Mz) = 0.118, 
/^soft = 2.5 GeV and Hhnvd = Mfac = 5 GeV, for various values of mini) and n = 4, 6, 8, 10, 20. The 
values for P*^ for rnijn) = 0.0 GeV are slightly different from the numbers shown in Ref. |^] as we 
used four-loop running and three-loop matching conditions at the five-four flavour threshold for the 
determination of strong coupling at the lower scales. In Ref. |^] two-loop running was employed. 



observations for the behaviour of the moments in the IS mass scheme made in the analysis of Ref. 
remain qualitatively true if the effects of light quark masses are included. Therefore we do not repeat 
any of the detailed discussions of Ref. |^] . 

In Table the values of at LO, NLO and NNLO have been displayed for M^^ = 4.7 GeV, 
ah'\Mz) = 0.118, /^soft = 2.5 GeV and //hard = /^fac = 5 GeV. for various values of m(m) and 
n = 4,6,8,10,20. 

Let us first discuss the results for n < 10. The numbers show that the light quark mass cor- 
rections are positive at NLO and negative at NNLO. (The observation that the NLO light quark mass 
corrections to the moments are positive is consistent with the NLO results for the charm quark mass 
effects in the perturbative bottom-antibottom IS wave function at the origin that were determined 
in Ref. [^4|.) The NNLO light quark mass corrections are actually larger than the NLO ones and 
lead to an overall negative shift in the moments at NNLO. For n = 4 the overall shift is around — 1% 
and for n = 10 around —5%, for m(m) ~ 1.5 GeV. This means that the bottom IS mass that is 
extracted from the T sum rule analysis receives a negative shift of the order of 15 MeV from the 
charm mass effects. This can be estimated from the overall mass dependence of the moments, see 
for instance Eqs. (|108| ) and (113). We will see in Sec. that this is indeed the case. Interestingly, 
for the determination of the bottom MS mass Mb(Mb), this means that the charm mass corrections 
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in the moments of the T sum rules and in the MS-IS mass relation are additive. The fact that the 
NNLO light quark corrections are larger than the NLO ones and lead to an overcompensation is not 
necessarily a point of concern, because the light quark mass corrections should be viewed as a part of 
the full corrections. Their behaviour is essentially not affected by that due to the small overall size of 
the light quark mass corrections. Nevertheless, one can and should ask why the behaviour of the light 
quark mass corrections in the moments is so much worse than in the bottom MS~1S mass relation 
discussed in Sec. ^. To judge the situation properly, however, one also has to take into account that 
the convergence of the massless corrections is also much worse in the moments than in the bottom 
MS-IS mass relation. (Compare the numbers shown in Table |5| with those in Table ^.) The difference 
in the behaviour of the perturbative series comes from the fact that the moments also depend on 
the square of the bottom-antibottom wave function at the origin (through their dependence on the 
bottom-antibottom decay and production rate). It is well known that the non-relativistic expansion 
of the bound state wave function at the origin is much worse than for the binding energy. A complete 
theoretical understanding of the physical origin of this behaviour, which could subsequently lead to an 
improvement of the situation, has not been achieved yet. However, it seems obvious that the problem 
arises from the infrared region]^. Paradoxically, the wave function at the origin is not affected by the 
question of the quark mass definition, which is generally argued to represent the dominant issue as 
far as infrared sensitivity is concerned; whether we use the bottom IS mass or the bottom pole mass 
does essentially not affect corrections to the wave function as shifts to the energy in the Schrodinger 
equation (^) leave the wave function unchanged. The behaviour of the light quark mass corrections 
in the moments therefore seems to be just a reflection of the present situation, in particular because 
the light quark mass corrections are, as mentioned several times earlier, very sensitive to moments 
smaller than the light quark mass. The present situation is that the moments of the T sum rules have 
a perturbative expansion that is much worse than the perturbative expansion of the bottom MS-IS 
mass relation, although general arguments based on global duality tell that the infrared sensitivity of 
the moments should be smaller , and although the gluon condensate contributions to the moments 
are negligibly small compared to the ones in the T(1S) mass. For the fits that we carry out in Sec. ^ 
to determine the bottom IS mass, we will disregard this problem as was done in previous T sum rule 
analyses. 

Let us now discuss the results for n = 20. We emphasize that moments Pn for n > 10 are not 
employed in our analysis. The light quark mass corrections for rn(rn) ~ 1.5 GeV at NLO and NNLO 
are scaringly large individually (of the order of 50% of the full corrections in each order), but lead 
to a relatively moderate overall correction at NNLO (of the order of 15%). We cannot exclude that, 
for very large values of n, the neglected double-insertion contributions at NNLO are very strongly 
enhanced for some reason, but we believe that this is very unlikely. The behaviour of the light quark 
mass corrections to the moments at very large values of n therefore seem to reflect the increased 
infrared sensitivity of the moments if n is considerably larger than 10. This is in agreement with the 
general considerations made by Poggio, Quinn and Weinberg |45| about the required minimal size of 
the energy interval in the integration that defines the moments in Eq. ( |106[ ). Thus the use of moments 
with n much larger than 10 could potentially lead to large systematic uncertainties and should be 
avoided. 



^ Voloshin ||55|] (see also Ref. [^) made the observation that the relative gluon condensate correction is about three 
times larger for the square of the wave functions at the origin than for the energy levels. 
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7 Determination of the Bottom Mass from the T Sum Rules 



In this section we carry out the numerical analysis for the determination of the bottom quark IS and 
MS masses using the T sum rules. In the analysis we put the main focus on the effects coming from 
the finite charm quark mass. In the previous sections the effects of the up, down and strange quark 
masses have have shown to be negligible. They are not included in the analysis. In Sec. 7.1 we will 
briefly review the statistical procedure that is used in the sum rule analysis, and in Sec. 7.2 we show 



the results for the bottom quark IS mass. In Sec. |7.3| the results for the bottom quark IS mass are 
confronted with the mass of the T(1S) and the implications for non-perturbative corrections in the 
T(1S) mass are examined. In Sec. the bottom MS mass is derived from the IS mass result. 



7.1 Statistical Procedure 

To obtain numerical results for the IS mass we use a statistical procedure based on the x^-function 

x'(Ml^af)(Mz),/isoft,/Xhard,/^fac) = E ( ^'f " ) (^~' W ( " ) , (161) 

{n},{m} 

where {n} represents the set of n's for which the fit is carried out and is the inverse covariance 
matrix describing the experimental errors and the correlation between the experimental moments. (See 
Ref. for a detailed description of the covariance matrix.) We emphasize that the renormalization 
scale dependence of the theoretical moments implemented in is exactly as displayed in the formulae 
in Sees. |6.2| and |6.3| and that no optimization procedure has been carried out. The covariance matrix 
contains the errors in the T electronic decay widths, the electromagnetic coupling aem! and the 
continuum cross section rcont- The small errors in the T masses are neglected. The correlations 
between the individual measurements of the electronic decay widths are estimated to be equal to 
the product of the respective systematic errors given in Table g In our two previous analyses in 
Refs. [0, m we have varied the correlation between this product and zero, and we found that the 
results are insensitive to this choice. In order to estimate the theoretical uncertainties in the mass 
extraction, the renormalization scales /Ugoft, /"hard and fifac are varied randomly in the ranges 

1.5 GeV < /isoft < 3.5 GeV, 
2.5 GeV < ^hard < 10 GeV, 

2.5 GeV < //fac < 10 GeV, (162) 

and the following sets of n's 

{n} = {4, 5, 6, 7} , {7, 8, 9, 10} , {4, 6, 8, 10} (163) 

are employed. Because the strong coupling cannot be extracted from the sum rules based on the 
function in Eq. ( |161| ) with relative uncertainties smaller than 10% (see the corresponding results of 
Refs. 0, I) we fit M^^, taking a^\Mz) as an inputf^. For each choice of the renormalization scales 

For this work we have also carried out a simultaneous fit of M^^ and ai^\Mz)- The results are not discussed in 
this work. As in the previous analyses in Refs. ^ we found a relative uncertainty of 10% in ai''\Mz)- The results 
for were found to be equivalent to those obtained when ai"^^ (Mz) is taken as an input. 
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/"soft, Mhardi A*fac in the ranges (|162| ), each set of n's and each given value of ai^\Mz), M^^ is obtained 
for which is minimal. 

We note that the dependence of the theoretical moments on variations of the renormalization 
scales, and in particular of //soft) is rather large, see Table 3 of Ref. This strong dependence is a 
reflection of the bad convergence properties of the moments caused by their dependence on the square 



of the bottom-antibottom wave functions at the origin, which we have already discussed in Sec. 6.4. 
It turns out that the x^-function in Eq. ( |161| ) has a smaller renormalization scale dependence than the 
individual moments, because it puts higher weight on linear combinations of the moments {P^ — P^) 
for which the overall dependence on /isoft in the wave functions partially cancels. This feature arises 
from the fact that the scale sensitivity of the squared wave functions for low radial bottom-antibottom 
excitations roughly corresponds to the size and correlation of the experimental uncertainties of the 
respective electronic T decay rates. Since the inverse covariance matrix puts higher weight on those 
linear combinations of the moments {P^ — P^) for which the experimental uncertainties coming 
from the electronic decay widths cancel, this consequently leads to partial cancellation of the scale 
dependence. This means that the x^-function essentially fits the relative size of the moments (as a 
function of n for a given value of /^soft) and not so much their overall normalization. In particular, 
this means that the resulting estimate for the theoretical uncertainty in the bottom IS mass based on 
the scale dependence of the x^-function in Eq. (161) is smaller than an estimate that is only based 



on a single moment. The estimate based on the x^-function is more realistic if the true theoretical 
uncertainties of the moments (and of the squared wave functions at the origin for the low radial 
excitations in particular), which are arising from truncating the perturbative series, are correlated 
according to their scale dependence. The stability of the estimate for the bottom IS mass based on 
the x^-fnnction in Eq. (161) under inclusion of higher order perturbative corrections (see Figs. ^ and 



|6| in Sec. [T^ ) seems to support that this is the case. 
7.2 Results for the Bottom IS Mass 

In Figs. ^ and |6|, the allowed range for the bottom IS mass is displayed as a function of the strong 
coupling based on the NLO and NNLO theoretical moments, respectively, for the charm MS masses 
fn{m) = 1.1, 1.3, 1.5 and 1.7 GeV. As a comparison, the results for m{fn) = 0.0 and 0.15 GeV are 
also shown. Each dot represents the bottom IS mass for which the function x^ is minimal for given 
(random) values of Og and the renormalization scales and a given set of n's. For the running of the 
strong coupling we employed the four-loop beta function and three-loop matching conditions at the 
five-four flavour threshold. The experimental and statistical errors that are indicated by the vertical 
lines correspond to 95% CL. They are below 15 MeV for all dots displayed in Figs. |5| and ^. We note 
that the NNLO results do not depend on whether the NNLO double-insertion charm mass corrections 
in the relation between the bottom pole and IS masses in Eq. ( |160D are taken into account or not. 

Figures |5| show that at NLO the inclusion of the finite charm mass does not at all affect the 
results for the bottom IS mass. For all choices of the MS charm quark mass between 1.1 and 1.7 GeV 
the allowed range for M^^ is practically identical to the analysis where the charm mass is neglected. 
Assuming 

a(5)(Mz) = 0.118 ± 0.003 (164) 



for the strong coupling, which is motivated by the compilations given in Refs. |57, 58], we arrive at 
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Figure 5: Results for the allowed range of M^^ for given values of a's\Mz) (and the corresponding 
values of af^ {2.b GeV)) at NLO for different choices of the MS charm quark mass. The dots represent 



points of minimal foi' a large number of random choices within the ranges (|162| ) and the sets (163), 
and randomly chosen values of the strong coupling. Experimental errors at 95% CL are displayed as 
vertical lines. It is illustrated how the allowed range for Mfjs is obtained if 0.115 < af\Mz) < 0.121 
is taken as an input. 



M^^ = 4.70 lb 0.06 GeV from the fits to the NLO theoretical moments. Although we assume a smaller 
uncertainty for as than in our previous analysis of Ref . |^] , where we used af^ (Mz) = 0.118 ± 0.004, 
the result for the bottom IS mass is equivalent to the one obtained at NLO in Ref. This is because 
the correlation of the allowed range for the IS mass to the value of the strong coupling is rather weak 
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Figure 6: Results for the allowed range of M^^ for given values of ah\Mz) (and the corresponding 
values of af^ (2.5 GeV)) at NNLO for different choices of the MS charm quark mass. The dots represent 



points of minimal for a large number of random choices within the ranges (|162| ) and the sets (163), 
and randomly chosen values of the strong coupling. Experimental errors at 95% CL are displayed as 
vertical lines. It is illustrated how the allowed range for Mfjs is obtained if 0.115 < af\Mz) < 0.121 
is taken as an input. 



(and because we round the results to 10 MeV precision). In contrast to the NLO results shown in 
Ref. Q, Fig. ^ exhibits an instability in the allowed IS mass range for input vales of ai'\Mz) larger 
than 0.125. This is caused by the use of the four- loop running for the strong coupling, which leads to 
larger values of a<j at lower scales than the two-loop running that has been used in Ref. |^] . For the 
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extraction of the bottom IS mass using Eq. ( 164 ) as an input, this instabihty is irrelevant. 

Figures ^ show the results for the allowed bottom IS mass range obtained from the NNLO 
theoretical moments. The results have a number of remarkable features. Most prominently, the already 
weak correlation of the allowed bottom IS mass range with the strong coupling visible in the NLO 
analysis has completely vanished in the relevant strong coupling range (Mz) ~ 0.118 if the effects 
of the charm quark mass are included. This leads to a variation of the bottom IS mass slightly smaller 
than the analysis with zero charm quark mass. We also find that the central value for the bottom IS 



mass is about 20 MeV lower if the charm quark mass is taken into account. Taking Eq. (164) as an 
input we arrive at 

Ml"" = 4.69 ± 0.03 GeV (165) 

for the bottom IS mass from the NNLO theoretical moments, taking into account the charm MS 
mass m(m) = 1.4 it 0.3 GeV. We emphasize that the uncertainty of 30 MeV in Eq. (|165| ) does not 
depend at all on the uncertainty of the strong coupling. It is also interesting that the allowed range for 
the bottom IS mass is stabilized for (Mz) larger than 0.125 compared to the NLO analysis. We 
note that the fact that the bottom IS mass gets a negative shift from the finite charm quark mass is 
expected, because a light quark mass gives negative corrections to the NNLO theoretical moments (see 
the discussion in Sec. |6.4| ). Neglecting the charm quark mass, the result reads M^^ = 4.71 ±0.04 GeV, 
where the variation is slightly larger than in Ref. owing to the use of the four-loop beta function 
for as- 

7.3 Comparison with the Mass of the T(1S) Meson 

It is a quite interesting feature of our NNLO result for the bottom IS mass, M^^ = 4.69 it 0.03 GeV, 
that it could be used to constrain the non-perturbative effects in the mass of the T(1S) meson, 
M(T(1S)) = 9460.37 lb 0.21 MeV. We remind the reader that, because we only used moments P„ with 
n < 10, the bottom-antibottom quark dynamics in the sum rules can be treated as perturbative (i.e. 
that the hierarchy "M^ ^ M^,v ^ M^,v'^ ^ Aqcd" is satisfied). Thus, non-perturbative corrections 
to the moments can be reliably determined in an expansion in gluonic and light quark condensates of 
increasing dimensions]^ For the moments P„ with n < 10 the dominant non-perturbative corrections 
from the dimension-4 gluon condensate corrections are below the per mille level (see Eq. ( [110| )). We 
therefore assume that non-perturbative effects in general can be safely neglected in our sum rule 



analysis. This means that our result for the bottom IS mass in Eq. (165) is not affected at all by 
non-perturbative effects. On the other hand, the bottom quark IS mass is just half of the perturbative 
contribution of the mass of the T(1S) meson, assuming that the latter is also a perturbative system. 

If we assume that the T(1S) meson can be treated as a perturbative system (at least as far as 
the calculation of its mass is concerned), non-perturbative corrections in M(T(1S)) can also be reliably 
determined in an expansion in the condensates. However, the corrections caused by the condensates 
are much larger than for the sum rules and cannot be neglected as in the sum rule analysis. A precise 
and accurate determination of the bottom IS mass from the T sum rules could therefore be used to 
seriously test the hypothesis that the T(1S) is perturbative and, eventually, to extract more precise 
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We assume for now that the expansion in terms of condensates is valid for the moments if n is not too large, and 



disregard the discussions in Sees. 6.1 and 6.4 
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values for the condensates. From the NNLO sum rule result for the bottom IS mass m Eq. ( |165D , we 
find 

M(T(1S))]"°"'"* = M(T(1S)) -2Mf = 80 ± 60MeV (166) 

for the non-perturbative contributions in the T(1S) mass. We emphasize again that the result in 
Eq. ( |166| ) can only be interpreted as "the non-perturbative contributions in the T(1S) mass" based on 
the hypothesis that the T(1S) is a perturbative system. It is interesting that Eq. (|166|) seems indeed 
consistent with the leading dimension-four gluon condensate contribution obtained by Voloshin and 



Leutwyler |5^, |5g| : 



-1 non-pert 1872 Mu 71^ TT r, 

*«^<'«»1 = 1275 (m:zw 

Equation ( p.67D obtains corrections from higher-dimension condensates as well as higher order per- 
turbative corrections to the Wilson coefficients multiplying each of the condensates. The subleading 



dimension-6 condensate contributions have been calculated in Ref. |59|. The unknown perturbative 
corrections induce a very large scale uncertainty of the dimension-4 condensate contribution shown 
in Eq. ( [L67[) . Assuming the standard literature range {ag ) = 0.05 it 0.03 GeV^ the leading gluon 
condensate correction can range from 10 up to 300 MeV (for 0.2 < < 0.5), where large numbers 
paradoxically correspond to smaller values of the strong coupling. As said before, this range is con- 



sistent with Eq. ( 166 ), but at the present stage it does not provide any meaningful determination of 
{ttg ), let alone a test of the hypothesis that the T(1S) is perturbative. The calculation of higher 



order perturbative corrections in Eq. (|167D would be very helpful to clarify the situation. 



7.4 Results for the Bottom MS Mass 

The heavy quark IS mass is, by design, a short-distance mass adapted to systems where the heavy 
quark is very close to its mass shell, or, in other word, where the heavy quark virtuality is much smaller 
than its mass, — Mq <C Mq. Apart from heavy-quark-antiquark systems the IS mass definition 
has also been applied to B mesons ||2^, where the bottom quark virtuality is of order Aqcd-^q- For 
high-energy or low-energy processes, however, where the heavy quark virtuality is of order or 
larger, the MS mass definition is the preferred one. In this section we determine the bottom MS mass 
from the NNLO result for the IS mass in Eq. ( |165| ) using formula (|82[). 

In Figs. the dependence of the bottom MS mass M^,{M^,) on the charm quark MS mass 
rn[rn), the strong coupling al (Mz), the renormalization scale fj, (of the strong coupling entering the 
relation between the bottom MS and IS masses) and the bottom IS mass M^^ is displayed. In each 
plot the respective fixed parameters have been chosen as as (Mz) = 0.118, M^^ = 4.69 GeV and 
IjL = 4.69 GeV, if not stated otherwise. Figure |^ shows Mb(Mb) as a function of the charm mass for 
fj. = 6.7 (dashed line), 4.69 (solid line) and 2.7 GeV (dotted line). Figures 0b,c,d show Mb(Mb) as a 
function of /u, and ai^^(Mz) for m(m) = 0.0 GeV (solid lines), 1.1 GeV (dotted lines), 1.3 GeV 
(dash-dotted lines), 1.5 GeV (dashed lines) and 1.7 GeV (long-dashed lines). It is conspicuous that for 
m(m) > 0.4 GeV and ^ > 2.5 GeV the dependence of Mb(Mb) on all four parameters is approximately 
linear. This welcome feature allows for the derivation of the following handy approximation formula: 



Mb(M, 



4.169 GeV - 0.01 (m(m) - 1.4 GeV ) + 0.925 ( - 4.69 GeV ) 
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Figure 7: The dependence of the bottom MS mass Mb(Mb), determined from the bottom IS mass, 
on the MS charm quark mass (a), the renormahzation scale ^ (b), the bottom IS mass M^^ (c) and the 



strong couphng ah {Mz) (d). In each plot the fixed parameters have been chosen as as (Mz) = 0.118, 
M^^ = 4.69 GeV and /i = 4.69 GeV, if not stated otherwise. In (a) the result is displayed for 
/X = 6.7 (dashed line), 4.69 (solid line) and 2.7 GeV (dotted line). In (b,c,d) the result is displayed for 
m(m) = 0.0 GeV (sohd lines), 1.1 GeV (dotted lines), 1.3 GeV (dash-dotted lines), 1.5 GeV (dashed 
lines) and 1.7 GeV (long-dashed lines). 



9.1 ( ) (Mz) - 0.118 ) GeV + 0.0057 ( fi - 4.69 ) GeV 



(168) 



For m{m) > 0.4 GeV and fi > 2.5 GeV the difference between this approximation formula and the 
exact result is less than 3 MeV. Figures ^and Eq. (168) show that at present the dominant uncertainties 



arise from the errors in the strong coupling and the bottom quark IS mass. The error in Mb (Mb) 
induced by M^^ is roughly equal to the error in M^^, and the one induced by ag roughly amounts to x 
times 10 MeV for an uncertainty of x times 0.001 in ai^^(M^). At the present stage M^^ and ai^\Mz) 
each contribute about 30 MeV to the uncertainty of Mb (Mb). The uncertainty coming from the scale 
variation and from the error in the charm mass corrections, which arises from the uncertainty in m(m) 
and the neglected order m^/Mf, corrections in the bottom pole-MS mass relation, is much smaller. 
Varying the renormahzation scale by 1 GeV shifts Mb (Mb) by about 5 MeV. This variation reflects the 
remaining perturbative uncertainty contained in the order a^ relation between the bottom MS and IS 
masses. We assign a perturbative uncertainty based on a 3 GeV variation of the renormahzation scale, 
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which gives 15 MeV. We estimate the error in the charm mass corrections as 5 MeV, which includes 
the sources of uncertainties in the charm mass corrections mentioned before and mass corrections 
coming from the other light quarks. Adding all uncertainties quadratically, we obtain 

M^(M^) = 4.17 ± 0.05 GeV, (169) 

where the central value has been determined for M^"" = = 4.69 GeV, af\Mz) = 0.118 and 
m(m) = 1.4 GeV. The result is rounded to units of 10 MeV. 

It is instructive to compare the result with the one obtained if the charm quark mass were 
neglected entirely: starting from 4.71 ±0.04 GeV for the bottom IS mass we would arrive at Mb (Mb) = 
4.20 GeV ± 30MeV[fe,] ± 40MeY[5M^''] ± 15MeV[/i] = 4.20 ± 0.05 GeV, where the central value is 
obtained for Mb'^ = ^ = 4.71 QeV and ai^^(Mz) = 0.118. Thus, after rounding to units of 10 MeV, 
the overall effect of the finite charm quark mass on the bottom MS mass is a shift of —30 MeV. 

Let us briefly comment on the results for Mb(Mb) that we gave in Refs. |^, 38|, where the 
bottom IS mass was also used as an intermediate step to determine Mb (Mb), but where charm mass 
effects were neglected. In Ref. @ we have obtained Mb(Mb) = 4.20 ± 0.06 GeV. The central value 
was obtained for Mb^ = A* = 4.71 GeV and af'\Mz) = 0.118, not using formula (|82|) for m{m) = 0, 
but solving Eqs. (^) and (|5^) numerically for a given value of Mb^. The order terms were not 



included because at that time the order corrections in the pole-MS mass relation |12, 35| where 



not yet known. The agreement of the central value with the one mentioned in the previous paragraph 
is a coincidence. The 60 MeV uncertainty was estimated by assuming a 40 MeV perturbative error, a 
40 MeV error caused by the strong coupling (assuming ai'\Mz) = 0.118 ± 0.004) and a 30 MeV error 
for Mb'^ The result Mb(Mb) = 4.21 ± 0.07 GeV in the proceedings quoted in Ref. |3| was obtained 
from formula (|8|) for m(m) = 0, assuming Mb'^ = 4.73 ± 0.05 GeV, /i = 4.73 GeV and af^Mz) = 
0.118 ± 0.004. The value for the bottom IS mass was determined from Mb'^ = 1/2(Mt(is) ± A^°|\"Py*) 
assuming Ax(i5) ~ 100 MeV for the non-perturbative contribution in the T(IS') mass. After rounding 
to units of 10 MeV, the error was 10 MeV larger than the one obtained from formula ( [8^ ) (for 



m{m) = 0) used in this work because in Ref. |38] the corresponding formula was based on the numerical 



order corrections to the pole-MS mass relation from Ref. |^5|. Equation (82), on the other hand, 
is based on the analytic result of Ref. [12|. 



8 Summary 

In this work we have examined, for the first time, the effects of light quark masses on bottom quark 
mass determinations based on sum rules for the T mesons at NNLO in the non-relativistic expansion. 
The effects coming from the charm quark mass are of particular interest, because one of the relevant 
physical scales governing the dynamics of the bottom-antibottom quark pair described in the sum 
rules is the inverse Bohr radius MbO^, which is as large as the charm quark mass. Thus, a naive 
expansion in the charm quark mass is a priori not possible in the calculations for the sum rules and 
the extraction of the bottom quark masses. Based on a NNLO analysis our result is that the charm 
mass effects amount to a shift of —20 MeV in the bottom IS mass, the mass parameter that is extracted 
directly from the sum rules. Our final result for the bottom IS mass reads 

Mb's = 4.69 ± 0.03 GeV , (170) 
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where the error arises from renormahzation scale variations; the correlation to the value of the strong 
coupling is negligible. The bottom IS mass is a short-distance mass, designed to be used in systems 
where the bottom quark virtuality is small with respect to its mass. Except for non-relativistic 
bottom-antibottom quark systems, it can also be applied directly to B meson decays. For inclusive 
semileptonic rates, it leads to very well behaved perturbative expansions |2^. Starting from the 
result for the bottom IS mass in Eq. ( |170D , we have examined the light quark mass effects in the 
determination of the bottom MS mass. Here, after rounding to 10 MeV precision, we find another 
shift of —10 MeV coming from the finite charm quark mass. Our final result for the bottom MS mass 
reads 

M^(Mb) = 4.17 ± 0.05 GeV, (171) 

where the 50 MeV error includes the error in Ml"", in the strong coupling {af\Mz) = 0.118 ± 0.003), 
the perturbative error estimated from variations of the renormahzation scale /x and the uncertainties 
in the charm mass corrections. Equation ( |171| ) represents a determination of the bottom MS mass at 
order a^. A handy approximation formula that shows the dependence of My^{M^) on the previously 
mentioned parameters can be found in Eq. (|168| ). We note that the charm quark mass effects in 
the bottom IS mass obtained from the sum rule analysis and those obtained in the relation between 
the bottom MS and IS masses have the same sign and lead to an overall shift of about —30 MeV 
in M^{M^). We emphasize that, for all practical applications, the values for the bottom IS and MS 
masses should, in contrast to the pole mass, not be considered as order-dependent numbers, since they 
do not contain an ambiguity of order Aqcd- Therefore, the results shown in Eqs. ( |170| ) and ( |171| ) can 
be used as an input for calculations at any order of perturbation theory. 

The order of magnitude of the charm mass corrections in the bottom MS mass can be antic- 
ipated from general considerations because the fact that one cannot expand in mcharm/CAfbOs) in the 
bottom IS-pole mass relation, whereas it is allowed to expand in m^harm in the bottom pole-MS mass 
relation, leads to are non-analytic correction oc a^m^t,^^^ that is vr^-enhanced and governed by the 
renormahzation scale rricharm rather than M^. The enhancement of the linear charm mass term can be 
understood from the fact that terms that a non-analytic in the square of the light quark masses are 
sensitive to small (infrared) momenta, in analogy to the non-analytic fictitious gluon mass terms that 
are often used in standard renormalon analyses. We have explicitly checked that the charm quark 
mass effects are under control in the bottom MS-IS mass relation, i.e. they do not lead to an un- 
controllable higher order uncertainty that would be relevant at present (or in the foreseeable future). 
This is because at higher orders of perturbation theory the charm quark (with finite mass) effectively 
decouples. It was one of the main motivations of this work to study the subtle interplay of the two 
mechanisms just mentioned at NNLO in the non-relativistic expansion. 

In this work we have determined the complete light quark mass corrections for the bottom IS 
mass obtained from the T sum rules at NNLO in the non-relativistic expansion, and for the bottom MS 
at order a^, with the following exceptions: in the moments for the sum rules, we neglected the NNLO 
double-insertion corrections from the NLO static potential at second order Rayleigh-Schrodinger per- 
turbation theory, since they are suppressed with respect to the single-insertion contributions coming 
from the NNLO static potential at first order Rayleigh-Schrodinger perturbation theory. For the 
charm mass corrections in the bottom pole-MS mass relation we have only taken into account the 
linear charm quark mass terms oc mcharm and have neglected higher order terms (?TT.(.harm/-^b)""icharm 
for n > 1. This is feasible because the bottom pole~MS mass relation does not involve the inverse 
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Bohr radius as a dynamical scale and because the subleading terms are less sensitive to small momenta 
and, therefore, not enhanced. We demonstrated that the linear charm mass approximation deviates 
from the full result for the charm mass corrections by only 10% for the cases where the full result is 
known. Finally, we have neglected the order light quark mass corrections oc a1{m^i,^,.^/AI^)'^ for 
n > 2 in the short-distance Wilson coefficient of the NRQCD bottom quark currents that describe 
bottom-antibottom production and annihilation in the non-relativistic regime. The neglected terms 
are small because the short-distance coefficient only contains contributions from momenta of order M^, 
which means that an expansion in the charm quark mass is justified. We estimate that all charm mass 
corrections that have been neglected in this work are an order of magnitude smaller than the ones that 
have been taken into account. We estimate the effects of the neglected terms and the uncertainties 
coming from the error in value of the charm quark mass to be at most at the level of 5 MeV. 

Note added: After this paper was completed we received Ref. where the N^LO large-/?o 
corrections in the heavy quark pole-lS mass relation where presented. The results agree with ours in 
App. g. 
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A NNLO Light Quark Mass Corrections to the S-Wave Heavy 
Quark— Antiquark Ground State Mass 

In this appendix we present details of the calculation of the light quark mass corrections to the ground 
state mass of the heavy-quark-antiquark J^*-" = ^Si bound state at NNLO in the non-relativistic 
expansion. The result is presented using the pole mass definition for the heavy and the light quark 
masses. We assume that the dynamics of the heavy-quark-antiquark pair is perturbative, i.e. that 
the scales Mq, MqV and MqV^ are assumed to be larger than Aqcd- We also assume that there are 
ni light quarks, of which one quark species has a finite pole mass m. We adopt the definition of the 
strong coupling where all n/ light quarks species contribute to the running. 

The formal expressions for the light quark mass corrections up to NNLO have been presented 
in Eqs. (|6l| ) and (]6^). In configuration space representation the expressions shown in Eqs. (61) and 
(H) read {x = \x\, y = \y\): 

-2M--A-°^„ = ld^x<p,six)SVZ^°ix)cP,six), (172) 
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■ 2 MP°'° A'~''~'^° 

(4 massive 



2 /d^y<^.s(x)<^K™(^)G,s(x,y)KZ^.L.(y)'A.s(y), (173) 



where K!L^°.iess 



^K!"^^:^!^): and dV^^^'-^ir) are given in Eqs. @, (0) and (|0|), respectively. The 
term (pis is the LO ground state wave function (see Eq. (|179D ) and 



Gisix,y) 



Ei — Eis 



y, 



(174) 



S-wavCji^l 

is the S-wave component of the configuration space LO Coulomb Green function Gc{x,y, E) for 
E = E^s = -2M^°'°A^°, where the ground state n = 1, ^Si energy-pole is subtracted. For the 
calculation of the NNLO light quark mass corrections only the sum over intermediate S-wave states 
has to be considered, because the light quark mass corrections to the static potential are spin- and 
angle- independent. An explicit expression for Gis{x,y) can be easily derived from the S-wave {I = 0) 
contribution of the configuration space Coulomb Green function representation by Voloshin ||5^ (a^ = 



G[='{x,y,E) 



M g Li__,{2kx)Li_,{2ky) 



in 



n=l 



[n 



1 



) n 



where 
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2 

fc2 = -ME 

and are the generalized Laguerre polynomials. The result reads 



(175) 

(176) 
(177) 



G'„s(x,y) 



lim 



G[=^ix,y,E) 



0ns(a;)(/'ns(y) 

EnS — E 



M7 / 7 



2 7rn 



exp 
27 



'27 



27 



'27 



2n^ 

00 



^(x + y)) lI_J—x) Ll_J—y) 

-^m-i(ir ^) -^m-i(ir y) 

(m — n)m 



for the subtraction of the n^Si energy pole, where 



vr 



E, 



-1/2^ 5/2^1 



MClal 



2^) Li_,(-^x 



nS 



4n2 



(178) 

(179) 
(180) 
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For n = 1 the first term in the curly brackets vanishes because L'^i{x) = 0. 

In the fohowing we present a number of intermediate results that lead to Eqs. (|6^ ) and 
All results are eventually expressed in terms of the function hi (i = 0, . . . , 7), which have been defined 
in Eqs. (|6|)-(|7D. 



Light Quark Mass Corrections at NLO 



The light quark mass corrections at NLO have been determined before in Ref. |£l[|. They consist of a 
single insertion of SV^^^ (Eq. (18)) in first order Rayleigh-Schrodinger perturbation theory: 



(^2 n2 



dx 



f{x) ] f ''^ \ 



11 

6" 



where 



(181) 



m 

7 



The function /iq can be easily calculated analytically, see Eq. (iGq). 



(182) 



Light Quark Mass Corrections at NNLO — Single Insertions 



For the light quark mass corrections at NNLO let us first consider the contributions coming from 
the single insertion of the NNLO potential at first order Rayleigh-Schrodinger perturbation theory, 
the first term on the RHS of Eq. ( |173| ). We already mentioned in the text after Eq. ( |30| ) that the 
light quark mass corrections to the NNLO static potential can be divided into three different parts, of 
which each vanishes individually for m ^ 0. In the following we present the results from these three 
contributions separately in terms of the functions hi: {m = e^^ m, fl = e"^^ ^), 



Cp as 



dx f{x) e 



-2mrx 



/5o(ln 



Ei(2 mxr) — Ei(— 2 mxr)^ — ai^ 



+ 3( \n[mr)+^- \ ("/^o ln(/ir) + ^"j + /3o ^ 



(^2 n2 

o Q 



3 / /47^\ \ 

- /io(a) In (-^j - flij - 3/3o/ii(o) 



3 . _ . 

/x2 

a\ ll\/„,/Ai\ „ ai\ 3 



(183) 



Cp as 



dx fix) + + (1 + 2x2) In "^^J 

\3 X2 V 2X \x+y/x^-l 



oo 

+ J dx fix) e-^™*^^ ln(4x2) - Ei(2 mxr)- Ei(-2 mxr)-^^ 
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5\2 7r2 



+ ( In(mr) + - 1 + 



o Q 



d^r<pUr) 



Cp as 



h2{a) + ho{a) ( 2 In + - 2h{a) + (ln[j) + f ] + ^ - 1 



57 / X61 13 

— ( ci r(0, 2 C2 m r) + di r(0, 2 m r) + ln(m ^) + ^ + Cs 



(184) 



o Q 



— — h — — ^ h ci ln(l + C2 a) + di ln(l + ^2 a, 



(185) 



Light Quark Mass Corrections at NNLO — Double Insertions 



For the light quark mass corrections at NNLO coming from double insertions of the NLO static poten- 
tial at second order Rayleigh-Schrodinger perturbation theory, there are two different contributions: 
the first contains two insertions of the NLO light quark mass corrections to the static potential, and 
the second contains one insertion of the NLO light quark mass corrections to the static potential and 
one insertion of the NLO static potential for massless light quarks. In detail, the two contributions 
read: 



^o(o) - ^ ho{a) + 2/13(0) - 2/14(0)^ - /15(a) + /i6(a) 



r2 ^2 



3 ho{a) ( In (^) + I) + 2 h{a) + 2 h,{a) ( In (0 + ^ 



1 ln2 

2 V2/ 



5 ^ /a\ 7r2 61 

— , , - In - - C3 H 

2 V2/ 6 V2/ ^ 6 72 



(186) 



2 /d^^ /d^y,^,s(x)<5Kr°(^)G,s(x,y)KriLe(y)</'is(y) 

/3o (^ho{a) (^3 In(^) -2) +2/i7(a) + 2/i3(a) (^l-ln(^)) 



^2 ^2 



+ ^ai(-^/.o(a) + /.3(a)-iln(|)-A 
Adding the two results we arrive at the expression displayed in Eq. 
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B Perturbative Heavy- Quark— Antiquark Ground State Mass at 
N^LO in the Large-/3o Approximation 

In this appendix we calculate the mass of the ground state of the heavy-quark-antiquark system at 
N^LO in the non-relativistic expansion (or order in the upsilon expansion) in the large-/3o approxi- 
mation, in terms of the heavy quark pole mass. It is assumed that the dynamics of the heavy quark 
pair is perturbative, i.e. that Mq » MqOs ^ Mqo'^ ^ Aqcd- All light quarks are treated as massless. 
Up to NNLO in the non-relativistic expansion, the large-/3o corrections can be easily extracted from 



the corrections shown in Eq. (53) (for zero light quark masses) by taking the highest power of ni and 



supplementing the result by the replacement n; — |/3o- 

Up to N'^LO the heavy-quark-antiquark ground state mass in the large-/3o approximation can 
be parametrized in the following way: 



m: 



n=l 

QQ 



2MX 



/3o 



2M, 



pole 



1 - eA^ 



A NLO 



, A NNLO 



^ A NNNLO 
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where (a^ = ai"' "*(/«)) 



C jp cig 



(189) 



A NLO 

^/3o 



A NNLO 



■ l2 + - L + - Cs + 



24 ^ 48 



In I 



(190) 

(191) 
(192) 



To calculate A^^^'^^'-' we first need the N^LO heavy-quark-antiquark static potential in the large-/3o 
approximation. In momentum space representation the large-/3o corrections are particularly simple 
because they only involve insertions of massless quark one-loop vacuum polarizations into the gluon 
line defining the LO static potential, supplemented by the replacement ni — > — |/3o5 



p2 



1 3 



+ 



In configuration space representation, the large-/3o potential reads 



(27r)3 



Vao(p) exp(ipr) 



In 



(193) 



(194) 
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where (/i = e'^^ /U, r 



M) 



ln(/ir) + - 





m ur H — H 

' 6/ 12 



ln(/ir) + ^') +!L ('ln(/ir) + ^) +2C3 



(195) 
(196) 

(197) 

(198) 



It is now straightforward to derive the formal result for A^™^*-" using Rayleigh-Schrodinger time- 
independent perturbation theory up to third order: 

2 MP°'° A^7^° = ( IS I y^f^^NNLo 1 15 ) + 2 ^ ( 15 I 



Eis — Ei 



rriT^iSnT^iS 



15 I V^';^" i;"^^^';' ( + 2M-'^ A 



NLO ^ I ) ( I T/NLO 



,00 



15), 



-2 1^3^ |d3y0,s(x)V>7^°(^)G,s(x,y)y^';^°(y)0,s(y) 
+ Jd'x jSy |d3zc/>,s(x)F^r(^)Gis(x,y)(F/,^°(y) + 2M--A^„^° 
xG,s(y,^)V>7°(;z)0is(^), 



where 



(pis{r) 



(r I IS) = (lS|r; 



"7 r 



(199) 



(200) 



(201) 



is the LO ground state Coulomb wave function and Gis{x,y) the S-wave component of the Coulomb 
Green function, where the ground state energy pole is subtracted. The explicit form of G is given in 
Eq. ( |178| ). The reader should note the appearance of the NLO binding energy in the third term on 
the RHS of Eq. ( [199| ) , which arises for the first time at third order Rayleigh-Schrodinger perturbation 
theory (see e.g. Ref. ||6^). For the calculation of the integrals in Eq. ( |199| ) the following formulae are 



helpful (ac = Cp a^"'\fi), x = \x\,y = \y\, z = |z|): 
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Using Eqs. (|202D - (|210D it is straightforward to derive the result for A^'^'~'^°- 
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where 



(212) 



C Order MS-IS Mass Relat ion in the Large- /^o Approximation 



In this appendix we determine the order (or order in the upsilon expansion) relation between 
the heavy quark MS and IS masses in the large-/3o approximation. We assume that there are ni hght 
quarks that are all massless. 

The large- /3o approximation of the heavy quark MS-pole mass relation is known to all orders 
of perturbation theory. The all order formula based on the resummation of massless quark one-loop 
vacuum polarization insertions into the one- loop gluon line reads [^H (a^ = ai"''*(MQ)) 



Mq(Mq)-M, 



pole 



E 
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1 
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(213) 



(214) 



(215) 



We note that in the large-/3o approximation the actual definition of the heavy quark mass as the 
renormalization scale in Eq. ( |213| ) is irrelevant, because this only affects subleading corrections in the 
large-/?o approximation. Evaluating Eq. ( pi 3] ) up to order we arrive at 



Mq(Mq 
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The large-/3o approximation for the MS-IS mass relation up to order (e^) using the upsilon expan- 
sion (with e = 1) then reads 



r(2) 



NNNLO 

/3o 



+ 



(221) 



where the A's are given in Eqs. ( |189 )-( |T91 ) and (211). We note that, in contrast to the complete 
result, the inversion of Eq. (|18^ ) is particularly simple in the large-/3o approximation because higher 
powers of low order contributions can be dropped. This is because they do not contain the highest 
power of ni in any order of perturbation theory. 



D Three-Loop Light Quark Mass Corrections to the Heavy Quark 
Pole— MS Mass Relation from Vacuum Polarization Insertions 

In this appendix we present analytic results for the functions fp''^ fp^^ and fpp (see Sec. which 
occur in the order and corrections in the relation between the heavy quark MS and the pole 
mass, and which come from insertions of massive quark one-loop vacuum polarizations into the gluon 
line (r = m/Mq): 
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The function P has been defined in Eq. (|T5|). The resuh for fp'^^ has aheady been obtained in Ref. |34] 



The results for fp^^ and fpp are new. 
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E Some Formulae for Moments with Massless Light Quarks 



In this appendix the expressions for the constants Wp'^'^ are given. They have been taken from Ref. 0]. 
They arise in the charm quark mass corrections to the moments of the T sum rules if the charm quark 
is included in the evolution of the strong coupling. The constants read {p = 1, 2, 3, . . .): 
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